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This is a review article of our work on hysteresis, avalanches, and criticality. We provide an 
extensive introduction to scaling and renormalization-group ideas, and discuss analytical and 
numerical results for size distributions, correlation functions, magnetization, avalanche durations 
and average avalanche shapes, and power spectra. We focus here on applications to magnetic 
Barkhausen noise, and briefly discuss non-magnetic systems with hysteresis and avalanches. 
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I. INTRODUCTION 

Our group has devoted several years to the study of 
the dynamics of a simple model, the random-field Ising 
model at zero temperature. We do clever simulations in- 
volving the dynamics of billions of domains (figure^), and 
we develop sophisticated analytical methods for extract- 
ing and explaining the properties of this model. We do 
this because we believe our model - despite its dramati- 
cally simplified nature - may well describe the properties 
of Barkhausen noise in many real physical systems. 

We present here the results of our simulations and anal- 
ysis, together with an explanation of why we believe they 
should (or could) describe real experiments. Our argu- 
ments for the applicability of our model are based on 
renormalization group and scaling theories, first devel- 
oped to study continuous phase transitions in equilibrium 
systems. To a large extent, these theories can be seen as 
the underlying reason why many if not most theories of 
nature apply to the real world, and (more specifically) 
why different magnets share common features in their 
dynamics despite having microscopically rather different 
morphologies and energetics. 

Why is the noise in magnets interesting? As one ap- 
plies an external field to most magnetic materials, the 
magnetization changes through the nucleation and mo- 
tion of domain walls. This motion is not smooth, and 
the resulting jumps in the net magnetization is termed 
Barkhausen noise. These jumps, corresponding to the 
reorganization (or avalanche) of a region of spins, usu- 
ally span many decades in size: there will be many small 
avalanches of spins, and fewer and fewer avalanches of 
larger and larger sizes (what we call crackling noise). In- 
deed, the probability D of having a jump or avalanche 




FIG. 1 Cross section of all avalanches in a billion-spin 
simulation at the critical disorder (JJ). The white background 
is the infinite, spanning avalanche, which would not exist 
above R c . 



of a given size S often decreases as a power law 1 in 
the size of the avalanche D(S) ~ S~ T for a significant 
range of sizes (figures |21 and [3J| . Simple behavior (like 
a power law) emerging out of complicated microscopic 
dynamics is a good sign that something needs explana- 
tion! Power laws in particular are the defining signature 
of continuous phase transitions, and a large fraction of 
the statistical physics community focuses on problems 
involving these power laws: a new class of systems to 
study is exciting. Finally, similar power laws for crackling 
noise emerge in many other systems (yj), from flux-line 
avalanches in superconductors, to crumpling paper (0), 
to earthquakes (the Gutenberg-Richter law). Magnetic 
Barkhausen noise forms a manageable experimental pro- 
totype of a whole class of behavior. 

Why should crackling noise in magnets be comprehen- 
sible? Very small jumps in the magnetization, corre- 
sponding to avalanches on scales comparable to the mi- 
crostructure, will certainly depend upon the details of 
the individual magnetic material. Large jumps, spanning 
the entire sample, will depend upon the sample geome- 
try. Indeed, these extremes in sizes are not regions in 
which simple power laws should be observed experimen- 
tally. But these extremes are separated by many powers 
of ten in a typical experiment. Just as the complex mi- 



1 In our model, this law is obtained only at the critical field and 
at the critical disorder (figure [3] integrating over field yields a 
different power law Dj nt (S) ~ S~ \ T + a P°) at the critical disorder, 
sec equation 1101 
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FIG. 2 Avalanche size distributions (Q). The figure 
shows the number D(S) of avalanches of size S measured in 
an experiment (Q) measuring near the center of the hystere- 
sis loop, our (nucleation) model near R c , and a related front 
propagation model. Notice the straight lines on a log-log plot 
indicate power-law behaviors D(S) ~ S~ T . 



croscopic properties of the molecules in a fluid only affect 
the viscosity and density on long length scales, so also one 
might expect the complicated microstructure in a mag- 
net could be subsumed into a few constants when dealing 
with events of sizes much larger than the microstructure. 
The smooth power-law behavior suggests exactly this: 
something simple is happening on intermediate scales, 
independent of either the microscopic or macroscopic de- 
tails. 

What should be explicable? How ambitious can we be 
in our expectations of a successful theory? As we shall 
explain in this review, we cannot expect to be able to 
tell when and where to expect an avalanche. (Our tools 
hence will likely not be effective at predicting large earth- 
quakes.) However, a successful theory should predict sta- 
tistical averages of almost any quantity that is dominated 
by events on large length and time scales, up to certain 
overall parameter-dependent scales (analogous to viscos- 
ity and density for fluids). For example, the probability 
D(S, T, L, W, R, H) in our theory of having an avalanche 
of size S, duration T, long axis L, short axis W, at dis- 
order R and external field H should be predicted by the 
theory, up to an overall size scale S s , time scale T s , field 
scale H s and offset H c . and disorder scale R s and offset 
R c . Indeed, we shall see that functions like these take 
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FIG. 3 Avalanche size distributions for our model (Q), 
integrated over the hysteresis loop, for a range of disorders 
above R c — 2.16: at R c one gets a power law (the straight 
dashed line) Di nt (S, R c ) oc s~^ T+al35 \ The inset shows a scal- 
ing collapse of the data (equation llOl . In the main figure, the 
thin lines show the scaling function prediction for the associ- 
ated curves. The true power-law behavior only emerges very 
close to R c , but the universal scaling function governs a range 
of disorders spanning a factor of two above R c . Notice also 
that we get six orders of magnitude of scaling at 5% above 
R c , without any self-organization. 



scaling forms 
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D s is a normalization fac- 



where r = ~ fl ° , h = — Tr 
tor, and r + avz is another universal critical exponent. 
We shall derive these scaling forms from self-similarity 
or scale invariance: the avalanche dynamics is inherently 
the same when observed (say) at microns and hundreds 
of microns. This scale invariance allows us to write func- 
tions of N variables as power laws times functions of 
N — 1 variables; hence, functions of one variable become 
pure power laws, as for D(S) above. What makes these 
theories predictive is that both the power law r and the 
entire multivariable function 2? mu iti are universal: they 
are independent of the microscopic model (within large 
classes of systems) . If your theory is in the same univer- 
sality class as the experiment, it will predict all critical 
exponents and scaling functions. 

There are three important qualifications about these 
predictions. 

1. Analytic corrections. The various scales S s , T s , 
R Cl H c , . . . will not be constants, but will depend 
smoothly and analytically on the control parame- 
ters in a given system. This is a serious issue, as 
the function we are describing (-D mu iti above) will 
typically itself be a smooth and analytic function of 
its parameters except at a phase transition where 



something qualitatively changes in the behavior of 
the system. Indeed, it is precisely at these phase 
transitions when events of all scales arise, where we 
expect our theory to apply. These analytic correc- 
tions usually become less and less important as we 
confine our attention to large, slow events near the 
transition. 



2. Universality Classes. The power of the theory 
rests in the prediction that quite different experi- 
mental systems, or a theoretical model and an ex- 
periment, may be in the same universality class. 
We shall see that universality classes are stud- 
ied by considering a space of possible systems: if 
two systems flow towards the same fixed point as 
one coarse-grains to larger scales, then their long 
length scale behaviors must agree. In many cases, 
and for magnetic noise in particular, there can be 
more than one candidate theory. Our model which 
incorporates short-range interactions and nucle- 
ation of new domains has competitors which al- 
low only interface depinning (|5|; |6|) and which in- 
corporate long-range fields in a mean-field fashion 
B EJ El Djf 0)- For many properties, the pre- 
dictions of the various theories are not strikingly 
different: careful experiments may be needed to 
distinguish between the rival theories if it is not 
obvious whether nucleation or long-range interac- 
tions are relevant in a given experiment. 

3. Dynamic criticality. The static, equal-time 
properties of these models have historically been 
less fussy than the time-dependent behavior. More 
specifically, there will often be several dynamic uni- 
versality classes for each static universality class, at 
least in the well-studied cases of equilibrium con- 
tinuous phase transitions. We will see clear evi- 
dence that our theories are not predicting the dy- 
namical behavior within the avalanches properly, 
but that the various theories are rather successful 
at describing the distributions of avalanches and 
other static quantities. 



II. MODELS 

Several variations of the zero temperature random field 
Ising model have been proposed to explain the power laws 
in Barkhausen noise. They are differentiated on the basis 
of the presence of long range forces, and the details of the 
dynamics. Our summary follows reference 1)131) . 

The model is composed of a large number of 'spins' 
Si = ±1 on a cubic lattice (square lattice in two dimen- 
sions, hypercubic for dimensions D > 3). These spins 
represent domains or small regions of the material. The 
model is subject to an external field H, which varies 
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FIG. 4 Hysteresis loop and subloops. The magnetization 
in our model (for R > R c ), as the external field H is ramped 
up and down. Our focus will primarily be on the upper, outer 
loop as the external field is ramped from — oo to oo. 



slowly in time, 2 usually ramping upward from H = — oo 
to H = oo. The value of a spin Si represents whether that 
domain is aligned (+1) or anti-aligned (—1) with the fi- 
nal H = +oo. The energy function, or Hamiltonian, for 
the models is 



H = - J^ JnnSiSj - ^2 Hs i - J! hiSi 
nn i i 

3 cos 2 (0j. 



(2) 



-Jj- Si - 2_^ Jdipole ' 
,; {hj} 



r?. 



S;Sj . 



Here J nn is the strength of the ferromagnetic nearest 
neighbor interactions, hi is a random field representing 
the effects of compositional and morphological disorder, 
Ji n f is the strength of an infinite range demagnetizing 
field, Jg) and 9ij is the angle between the positive spin 
direction and the difference vector between lattice posi- 
tions i and j, and Jdipoic is the strength of the dipolc- 
dipole interactions. The critical exponents of the power 
laws are independent of the particular choice of random 
field distributions p(hi) for a large variety of distribu- 
tions. We use a Gaussian distribution of random fields, 
with zero mean and standard deviation R. (When we 
refer to the strength of the disorder, we are referring to 
the width, R, of the random field distribution.) 3 



Except when otherwise mentioned, the field varies slowly enough 
that each avalanche finishes before the field changes. 
In this review we focus specifically on zero temperature mod- 
els; ther e are al so finite temperature studies espe cially for pure 
systems ll4t Il5l) and micromagnetic systems {UjjL int erface de- 
pinning models in the presence of temperature Jl7t Lffi and hys- 
teresis due to driven interfaces at finite temperature (19; 20; 21). 



Two different dynamics have been considered. The 
first is a front propagation dynamics in which a spin on 
the edge of an existing front flips as soon as it would 
decrease the energy to do so. introduced by Ji and Rob- 
bing (|5|). Spins with no flipped neighbors cannot flip even 
if it would be e nerg etically favorable. Second is the dy- 
namics we use (23), which includes domain nucleation. 
Any spin can flip when it becomes energetically favor- 
able to do so. In both cases, spins flip in shells — all spins 
which can flip at time t flip, then all of their newly flip- 
pable neighbors flip at time t+1, causing an avalanche. 
The number of spins flipped in each shell gives the time 
series of the avalanche (figure [SJ, whose irregular fluctua- 
tions and near halts are typical also of experimental time 
series. 




FIG. 5 Typical avalanche time series. il3tl Voltage (num- 
ber of domains flipped) pulse during a single large avalanche 
(arbitrary units). Notice how the avalanche almost stops sev- 
eral times: if the forcing were slightly smaller, this large 
avalanche would have broken up into two or three smaller 
ones. The fact that the forcing is just large enough to on 
average keep the avalanche growing is the cause of the self- 
similarity: on average a partial avalanche of size S will trigger 
one other on size 5*. 

Depending on which terms are included in the Hamil- 
tonian, the behavior appears to fall into three different 
universality classes. 

1. Front propagation model (5). The front prop- 
agation model has a critical field H c for a range 
of disorders R. As the field approaches H c , larger 
and larger avalanches arise. As one approaches H c 
from below these avalanches develop a power-law 



We are currently working on hy ster esis at finite temperatures 
with randomness and nucleation 1 22). 



distribution in sizes with a cutoff which diverges at 
H c . Above H c the entire front is depinned, and the 
front moves forward with an inhomogeneous veloc- 
ity which begins jerky on all length and time scales 
(echoing the avalanches below H c ) but which be- 
comes smoother and more uniform at large H. 
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FIG. 6 Infinite avalanche, jump in M(H) for 71 < R c . 

2. Our model (nucleation allowed, Figure l4|l . 

For our model where domain nucleation is allowed 
and only nearest neighbor interactions are included, 
there is a continuous transition at a critical point 
with disorder R c and external field H c . At very 
large disorder, the couplings between spins J nn are 
a small perturbation and each spin s» flips roughly 
when the external field H(t) is equal and opposite 
to the local random field hi- For all disorders above 
R c this qualitative behavior persists: all avalanches 
are finite, dying out after a number of spin flips that 
does not depend on the number of spins N, yielding 
a smooth hysteresis loop (figure EJ). For very small 
disorder, an early spin will trigger almost the entire 
system to flip over. This too extends qualitatively 
to all disorders below R c , where finite fraction of 
the spins in the system (even as the system size 
becomes large) flip in a single event, the infinite or 
spanning avalanche. The critical disorder R c is de- 
fined as the disorder at which the infinite avalanche 
first arises, and where a jump in the magnetization 
per spin m(H) first arises in the limit of an infi- 
nite system. 4 Figure 2] shows the avalanches in a 
cross-section of our system at the critical disorder: 
the white background is the spanning, or infinite 



avalanche. 5 One does not need to be exactly at 
R c to observe power laws in our model: one gets 
a large region of power-law scaling rather far from 
the critical point (figure 03 and reference iQ)). 

3. Infinite Range and Dipolar Fields. For dis- 
orders below R c , or when domain nucleation is 
not allowed, the addition of an infinite-range de- 
magnetizing field 6 (@) can self-organize the system 
to a different critical behavior. (Self-organization 
means that the system naturally sits at a criti- 
cal point, without having to tune any parameters.) 
The infinite-ranged interaction is sometimes intro- 
duced to mimic the effects of the boundaries of ma- 
terials with dipolar, or other long-ranged interac- 
tions q lQf) : these interactions also self-organize the 
model to the critical point {j|). Our use of the term 
"infinite-range models" perhaps obscures the clear 
physical origin of this universality class of models. 

Zapperi et al. (flOf) argue that the addition of dipole- 
dipole interactions to the model lowers the upper 
critical dimension to three and produces mean- 
field exponents in three dimensions. Since large 
mean-field simulations are much easier than large 
simulations with dipole-dipole interactions, we will 
give results from mean-field simulations in this pa- 
per. Dipolar interactions without an infinite-range 
term were explored by Magni <J24l) in two dimen- 
sions, who found labyrinthine patterns and hystere- 
sis loops similar to those seen in garnet films. 

The infinite-range model is apparently also equiv- 
alent to the rather successful single-degree of free- 
dom ABBM models Gl El El • 



Unless we specify otherwise, we will focus on our 
model, allowing nucleation with parameters J nn , H, and 
hi but without long-range interactions (J; n f = Jdipoic = 
0). 



III. THE RENORMALIZATION GROUP AND SCALING 

To study crackling noise, we use renormalization-group 
tools developed in the study of continuous phase transi- 
tions. 7 The word renormalization has roots in the study 



4 Each avalanche produces a jump in the magnetization AM, but 
all finite avalanches produce jumps in the magnetization per spin 
=¥■ which vanish as the number of spins N diverges. 



It may seem as if the spanning avalanche is occupying quite 
a large fraction of the volume. As the system size increases, 
the "holes" formed by the finite avalanches (and other spanning 
avalanches) will gradually become larger and larger fractions of 
the total, and the spanning avalanche at R c should grow with 
the number of spins TV to a power below one. 
This is not the same as a mean-field model, because the model 
contains nearest neighbor interactions along with the infinite 
range interactions. 

This section and the next follow closely the presentation in ref- 
erence lj}). References to the broader literature may be found 
there. 



of quantum electrodynamics, where the effective charge 
changes in size (norm) as a function of length scale. The 
word group refers to the family of coarse-graining opera- 
tions basic to the method: the group product is compo- 
sition (coarsening repeatedly). The name is unfortunate, 
however, as the basic coarse-graining operation does not 
have an inverse, and thus the renormalization group does 
not have the mathematical structure of a group. 
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Infinite 
Avalanche 



FIG. 7 Renormalization— group flows ((J). The 

renormalization-group uses coarse-graining to longer length 
scales to produce a mapping from the space of physical sys- 
tems into itself. Consider the space of all possible systems 
exhibiting magnetic hysteresis (including, in our imagination, 
both real models and experimental systems) . Each model can 
be coarse-grained, removing some fraction of the microscopic 
domains and introducing more complex dynamical rules so 
that the remaining domains still flip over at the same exter- 
nal fields. This defines a mapping of our space of models into 
itself. A fixed point S* in this space will be self-similar: be- 
cause it maps into itself upon coarse-graining, it must have 
the same behavior on different length scales. Points that flow 
into S* under coarse-graining share this self-similar behav- 
ior on sufficiently long length scales: they all share the same 
universality class. 

The renormalization group studies the way the space 
of all physical systems maps into itself under coarse- 
graining (see figure ■ The coarse-graining operation 
shrinks the system, and removes degrees of freedom on 
short length scales. Under coarse-graining, we often find 
a fixed point S*: many different models flow into the 
fixed point and hence share long-wavelength properties. 
To get a schematic view of coarse-graining, look at fig- 
ures 0and OD the 1000 3 cross section looks (statistically) 
like the 100 3 section if you blur your eyes by a factor of 
10. 

Much of the mathematical complexity of this field in- 
volves finding analytical tools for computing the flow di- 
agram in figure Using methods developed to study 
thermodynamical phase transitions and the depinning of 



charge-density waves, we can calculate for our model the 
flows for systems in dimensions close to six (the so-called 
e expansion, where e — 6 — D, D being the dimension of 
the system) . Interpolating between dimensions may seem 
a surprising thing to do. In our system it gives rather 
good predictions even in three dimensions (i.e., D = 3), 
but it's hard work, and we won't discuss it here. Nor will 
we discuss real-space renormalization-group methods or 
series expansion methods. We focus on the relatively 
simple task of using the renormalization group to justify 
and explain the universality, self-similarity, and scaling 
observed in nature. 
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FIG. 8 Cross section of all avalanches in a million-spin 
simulation at the critical disorder (.1). Compare with the 
much larger simulation in figure 



Consider the "system space" for disordered magnets. 
There is a separate dimension in system space for each 
possible parameter in a theoretical model (disorder, cou- 
pling, next-neighbor coupling, dipolar fields etc.) or in an 
experiment (temperature, annealing time, chemical com- 
position ...). Coarse-graining, however one implements it, 
gives a mapping from system space into itself: shrinking 
the system and ignoring the shortest length scales yields a 
new physical system with identical long-distance physics, 
but with different (renormalized) values of the param- 
eters. We've abstracted the problem of understanding 
crackling noise in magnets into understanding a dynami- 
cal system acting on the space of all dynamical systems. 

Figure represents a two-dimensional cross section of 
this infinite-dimensional system space. We've chosen the 
cross section to include our model: as we vary the dis- 
order i?, our model sweeps out a straight line in system 
space. The cross section also includes a fixed point S*, 
which maps into itself under coarse-graining. The sys- 
tem S* looks the same on all length and time scales, 



because it coarse-grains into itself. We can picture the 
cross section of figure \7\ either as a plane in system space 
(in which case the arrows and flows depict projections, 
since in general the real flows will point somewhat out 
of the plane), or as the curved manifold swept out by 
our one-parameter model as we coarse grain (in which 
case the flows above our model and below the horizontal 
curved line in figure {7\ should be ignored). 

The flow near S* has one unstable direction, leading 
outward along the horizontal curve (the unstable man- 
ifold). In system space, there is a surface of points C 
which flow into S* under coarse-graining. Because S* 
has only one unstable direction, C divides system space 
into two phases. To the left of C, the systems will have 
one large, system-spanning avalanche (a snapping noise). 
To the right of C, all avalanches are finite and under 
coarse-graining they all become small (popping noise). 
Our model, as it crosses C at the value R c , goes through 
a phase transition. 

Our model at R c is not self-similar on the shortest 
length scales (where the square lattice of domains still is 
important), but because it flows into S* as we coarse- 
grain we deduce that it is self-similar on long length 
scales. Some phase transitions, like ice melting into wa- 
ter, are abrupt and don't exhibit self-similarity. Con- 
tinuous phase transitions like ours almost always have 
self-similar fluctuations on long length scales. Also, we 
must note that our model at R c will have the same self- 
similar structure as S* does. Indeed, any experimental 
or theoretical model lying on the critical surface C will 
share the same long-wavelength critical behavior. This is 
the fundamental explanation for universality. 

The flows in system space can vary from one class of 
problems to another: the system space for front propa- 
gation models (figure EK) will have a different flow, and 
its fixed point will have different scaling behavior (yield- 
ing a different universality class). In some cases, a fixed 
point will attract all the systems in its vicinity (figure^J . 
Usually at such attracting fixed points the fluctuations 
become unimportant at long length scales: the Navier- 
Stokes equation for fluids described earlier can be viewed 
as a stable fixed point. The coarse-graining process, av- 
eraging over many degrees of freedom, naturally smooths 
out fluctuations, if they aren't amplified near a critical 
point by the unstable direction. Fluctuations can remain 
important when a system has random noise in a con- 
served property, so that fluctuations can only die away 
by diffusion: in these cases, the whole phase will have 
self-similar fluctuations, leading to generic scale invari- 
ance. 

Sometimes, even when the system space has an un- 
stable direction like in figure Q the observed behavior 
always has avalanches of all scales. This can occur sim- 
ply because the physical system averages over a range 
of model parameters (i.e., averaging over a range of R 
including R c in figured). For example, this can occur 
by the sweeping of a parameter slowly in time, or vary- 
ing it gradually in space either deliberately or through 




FIG. 9 Flows for a front —propagation model (f). The 
front propagation model has a critical field H c at which the 
front changes from a pinned to sliding state, (a) Coarse- 
graining defines a flow on the space of depinning models. The 
fixed point S*/ p will have a different local flow field from other 
renormalization-group fixed points, yielding its own univer- 
sality class of critical exponents and scaling functions. The 
critical manifold C, consisting of models which flow into S*/ p , 
separates the stuck fronts from those which slide forward with 
an average velocity v(H). (b) The velocity varies with the 
external force as a power law v(H) ~ (H — H c ) ■ Clever 
experiments, or long-range fields, can act to control not the 
external field, but the net magnetization: changing the mag- 
netization slowly sets v ~ 0, thus self-tuning H « H c - This 
is one example of self-organized criticality. 



large-scale inhomogeneities. 

One can also have self-organized criticality, where the 
system is controlled so that it naturally sits on the crit- 
ical surface. Self-organization to the critical point can 
occur via many mechanisms. In the front propagation 
models with long-range fields (figure Et>), upon an incre- 
ment of external field H the front will move forward until 
the long-range demagnetization force balances the incre- 
ment: the net external field given by the sum organizes 
itself to precisely the critical depinning field for the front. 



IV. CRITICAL EXPONENTS, SCALING FUNCTIONS 

The renormalization group is the theoretical under- 
pinning for understanding why universality and self- 
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FIG. 10 Attracting fixed point (1). Often there will be 
fixed points that attract in all directions. These fixed points 
describe phases rather than phase transitions. Most phases 
are rather simple, with fluctuations that die away on long 
length scales. When fluctuations remain important, they will 
exhibit self-similarity and power laws called generic scale in- 
variance. 



similarity occur. 8 Once we grant that different systems 
should sometimes share long-distance properties, though, 
we can quite easily derive some powerful predictions. 

To take a tangible example, lets consider the relation 
between the duration of an avalanche and its size. If we 
look at all avalanches of a certain duration T in an experi- 
ment, they will have a distribution of sizes S around some 
average (S) oxp (T). If we look at a theoretical model, it 
will have a corresponding average size (S}th(T). If our 
model describes the experiment, these functions must be 
essentially the same at large S and large T. We must al- 
low for the fact that the experimental units of time and 
size will be different from the ones in our model: the best 
we can hope for is that 



(SWT) = A(S) th (T/B), 



(3) 



for some rescaling factors A and B. 

Now, instead of comparing to experiment, we can com- 
pare our model to itself on a slightly larger time scale. If 
the time scale is expanded by a small factor B = 1/(1— 5), 
then the rescaling of the size will also be small, say 1 + a. 
Hence 



(S)(T) = (l+a6)(S)((l-6)T). 



(4) 



Making 6 very small yields the simple relation a (S) = 
Td(S)/dT, which can be solved to give the power law re- 
lation (S)(T) = SoT a . The exponent a is called a critical 
exponent, and is a universal prediction of a given theory. 



(That means that if the theory correctly describes an ex- 
periment, the critical exponents will agree.) In our work, 
we write the exponent a relating time to size in terms of 
three other critical exponents, a — 1/avz. 
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FIG. 11 Universal Critical Exponents vs. Experi- 
ment. (1) Different experiments on crackling noise in mag- 
nets measure different combinations of the universal critical 
exponents. Here we compare experimental measurements (see 
table I of reference Q)) to the theoretical predictions for three 
models: our model, the front-propagation model and the dipo- 
lar mean-field theory. Power laws giving the probability of 
getting an avalanche of a given size, duration, or energy at 
the critical point are shown; also shown is the critical expo- 
nent giving the power as a function of frequency (Ua) (due to 
the internal structure of the avalanches, figure|^J . In each pair 
of columns, the first column includes only avalanches at exter- 
nal fields H near H c where the largest avalanches occur, and 
the second column (when it exists) includes all avalanches. 
The various combinations of the basic critical exponents can 
be derived from exponent equality calculations similar to the 
one discussed in the text. Many of the experiments were done 
years before the theories were developed: many did not report 
error bars. All three theories do well (especially considering 
the possible systematic errors in fitting power laws to the ex- 
perimental measurements: see in figure|^|how a bit away from 
R c the effective slopes change). (A more systematic and criti- 
cal review of the exponent measurements in the literature may 
be found in Table I of reference l|ll() . One may also find in 
section V. A of their paper a discussion of the relation between 
the experiments and the theoretical universality classes.) 

There are several basic critical exponents, which arise 
in many different combinations depending on the phys- 
ical property being studied. We've seen that at R c and 
near H c the probability of having an avalanche of size S 
goes as S~ T . 9 The cutoff in the avalanche size distribu- 
tion in figure |3 gets larger as one approaches the critical 



This section also follows closely the presentation in reference y|) . 



9 We will see that when the entire hysteresis loop is considered this 
changes to S- < - T + a ^ s '>. 



disorder as (R — Rc)~ a (figure |3J). The typical length of 
the largest avalanche goes as (R — R c )~ v - The jump in 
the magnetization (ngurelHJ goes as (R — R C )P, and at R c 
the magnetization (M - M c ) ~ (H - H c ) l / S . 10 Finally, 
we need to know how time rescales: the duration of an 
avalanche of spatial extent L typically will go as L z . 

Any physical property that shows singular behavior at 
the critical point will have a critical exponent that can 
be written in terms of these basic ones ( figure ITU) . 11 

To specialists in critical phenomena, these exponents 
are central; whole conversations will seem to rotate 
around various combinations of Greek letters. Critical 
exponents are one of the relatively easy things to calcu- 
late from the various analytic approaches, and so have 
attracted the most attention. They are derived from the 
eigenvalues of the linearized flows about the fixed point 
S* in figure [7| Figure IT2l shows our numerical estimates 
for several critical exponents in our model in various spa- 
tial dimensions, together with our 6 — e expansions for 
them. Of course the key challenge is not to get analytical 
work to agree with numerics: it's to get theory to agree 
with experiment. Figure ITU shows that our model does 
rather well in describing a wide variety of experiments, 
but that the two rival models (with different flows around 
their fixed points) also fit. 

Critical exponents are not the be-all and end-all: many 
other scaling predictions, explaining wide varieties of be- 
havior, are quite easy to extract from numerical simu- 
lations. Universality extends even to those long length 
scale properties for which one cannot write formulas. 
Perhaps the most important of these other predictions 
are the universal scaling functions. For example, lets con- 
sider the time history of the avalanches, V(i), denoting 
the number of domains flipping per unit time. (We call it 
V because it's usually measured as a voltage in a pickup 
coil.) Each avalanche has large fluctuations, but one can 
average over many avalanches to get a typical shape. Fig- 
ures ED and [TI| show averages over all avalanches of fixed 
duration T. Lets call this (V)(T,t). Universality again 
suggests that this average should be the same for experi- 
ment and a successful theory, apart from an overall shift 
in time and voltage scales: 
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FIG. 12 Universal Critical Exponents in Various Spa- 
tial Dimensions. Q) We test our e-expansion predic- 
tions 11271) by measuring the various critical exponents nu- 
merically in up to five spatial dimensions (y; 1281) . The var- 
ious exponents are described in the text. All of the ex- 
ponents are calculated only to linear order in e, except for 
the correlation length exponent v, where we use results from 
other models. The agreement even in three dimensions is re- 
markably good, considering that we're expanding in 6 — D 
where D = 3! We should note that perturbing in dimen- 
sion for our system is not only complicated, but also contro- 
versial. Our ex pan sion uses the Martin Siggia Rose formal- 
ism dH ||3 mill |H Hi |H HI) to describe a determinis- 
tic dynamical system without any temperature fluctuations, 
while the calculation for the equilibrium model involves tem- 
perature fluctuations and no history dependence at all. Even 
so, we have shown that our 6 — e expansion for the critical ex- 
ponents of the zero-temperature, nonequilibrium model maps 
to all orders in e onto that of the equilibrium model. The 
6 — e expansion for the equilibrium model, however, has been 
controversial for decades l|37l). Recently the thermal 6 — e 
expansion of the equilibrium model has been called into ques- 
tion (37; 38; 39) not just due to non perturbative corrections, 
as was previously assumed, but due to previously neglected 
higher order terms in the expansion. The implications of this 
controversy for our non-equilibrium renormalization group 
treatments is not yet known. 



(U) oxp (T, t) = A{V) th (T/B, t/B). 



(5) 



Comparing our model to itself with a shifted time scale 
becomes simple if we change variables: let v(T,t/T) = 
(V)(T,t), so v(T,t/T) = Av(T/B,t/T). Here t/T is a 
particularly simple example of a scaling variable. Now, 
if we rescale time by a small factor B — 1/(1 — S), we 
have v(T,t/T) = (1 + b)v(t/T, (1 - S)T). Again, making 
5 small we find bv = Tdv/dT, with solution v = VoT b . 



10 Don't confuse the small change in scale <5 earlier with the critical 
exponent S here. 

11 Indeed, there are relations even between these exponents: see 
section lVll 



However, the integration constant vq will now depend on 
t/T, vQ = V(t/T), so we arrive at the scaling form 



(V)(t,T) = T b V(t/T), 



(0) 



where the entire scaling function V is a universal predic- 
tion of the theory. 

Figures ITS1 and ITU show the universal scaling functions 
V for two models and three experiments. For our model, 
we've drawn what are called scaling collapses, a simple 
but powerful way to both check that were in the scaling 
regime, and to measure the universal scaling function. 
Using the form of the scaling equation Eq. [2] , we simply 
plot T b (V)(t,T) versus t/T, for a series of long times T. 
All the plots fall onto the same curve. This tells us that 
our avalanches are large enough to be self-similar. (If in 
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FIG. 13 Comparison of experimental average pulse 
shapes for fixed pulse duration, as measured by three dif- 
ferent groups (Qill3iIl2iE3)- Our theories don't only predict 
power laws: they should describe all behavior on long length 
and time scales (at least in a statistical sense). In particular, 
by fixing parameters one can predict what are called scaling 
functions. If we average the voltage as a function of time over 
all avalanches of a fixed duration, we get an average shape. In 
these three experiments, they find that this shape is the same 
for different durations: by rescaling the time by the dura- 
tion T and the voltage by the maximum of the average curve, 
the curves within a given experiment collapse onto one. The 
three experiments, however, do not all collapse onto the same 
curve. This could mean that they are in different universality 
classes, or that we don't understand this dynamical scaling 
completely. 



your scaling collapse the corresponding plots do not all 
look alike, then any power laws you have measured are 
probably accidental.) It also provides us with a numerical 
evaluation of the scaling function V. Note that we use 
1/avz — 1 for the critical exponent b. This is an example 
of an exponent equality: easily derived from the fact that 



(S)(T) = / (V)(t, T)dt= / T b V(t/T)dt ~ T 



nb+i 



(7) 



and the scaling relation (S)(T) ~ T- X l avz . 

Notice that the two models and the three experiments 
have quite different shapes for V. How do we react to 
this? Our models arc falsified if any of the predictions 
are shown to be wrong asymptotically on long length 
and time scales: hence our theory is either wrong (in- 
applicable to these systems) or somehow at least incom- 



FIG. 14 Comparison of theoretical average pulse shape 
scaling functions for our nucleated model and the front 
propagation model (Q)- The front propagation models have 
1/avz — 1.72 ± 0.03 in this collapse; our nucleation model 
has \javz = 1.75 ± 0.03 (in principle, there is no reason to 
believe these two should agree). The inset shows the two 
curves rescaled to the same height (the overall height is a non- 
universal feature): they are quantitatively different, but far 
more similar to one another than either is to the experimen- 
tal curves in figure [TBI The mean-field model apparently has 
a scaling function which is a perfect inverted parabola (141ft . 
The ABBM model, which is also mean-field like, interestingly 
has a different shape, that of one lobe of a sinusoid <I421 eq. 
102). 



plete. Incorporating insights from careful experiments to 
refine the theoretical models has historically been crucial 
in the broad field of critical phenomena. The message 
we emphasize here is that scaling functions can provide 
a sharper tool for discriminating between different uni- 
versality classes than critical exponents. 

Broadly speaking, most common properties that in- 
volve large length and time scales have scaling forms: 
using self-similarity, one can write functions of N vari- 
ables in terms of scaling functions of N — 1 variables: 
F(x,y,z) = z~ a F(x/z l3 ,y/z~ f ). In the inset to figure |3 
we show the scaling collapse for the integrated avalanche 
size distribution (derived as equation II Oj l: 

D int (S,R) = S-( T +^V int ((R~R c )/S-°). (8) 

This example illustrates that scaling works not only at R c 
but also near R c ; the unstable manifold in figure \J\ gov- 
erns the behavior for systems near the critical manifold 
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c. 



V. MEASURING EXPONENTS AND SCALING 
FUNCTIONS 

The simulations provide a rich variety of measured 
quantities, each with a characteristic universal scaling 
form and associated critical exponents. In this section 
we focus on a few key ones of particular experimental 
significance. 12 

We'll discuss the following properties obtained from 
the simulation 

• the magnetization M(H, R) as a function of the 
external field H. 

• the avalanche size distribution integrated over the 
field H, D int (S,R). 

• the avalanche correlation function integrated over 
the field H, Gi nt {x, R). 

• the distribution of avalanche durations D\ m ' (S, t) 
as a function of the avalanche size S, at R = R Cl 
integrated over the field H . 

• the energy spectrum E(lj) of the Barkhausen noise 
for various models at criticality. 



a 
a 

N 

a 

bo 




0.0 



S -1-0 



-2.0 



1.0 1.5 

Magnetic Field (H/J) 



2.0 




(h+Br)/lrl 



A. Magnetization Curves 

Unfortunately the most obvious measured quantity in 
our simulations, the magnetization curve M(H), is the 
one which collapses least well in our simulations. We 
start with it nonetheless. 

The top figure 1151 shows the magnetization curves ob- 
tained from our simulation in 3 dimensions for several 
values of the disorder R. As the disorder R is decreased, 
a discontinuity or jump in the magnetization curve ap- 
pears where a single avalanche occupies a large fraction of 
the total system. In the thermodynamic limit this would 
be the infinite avalanche: the largest disorder at which it 
occurs is the critical disorder R c . For finite size systems, 
like the ones we use in our simulation, we observe an 
avalanche which spans the system at a higher disorder, 
which gradually approaches R c as the system size grows. 

Figure LSI shows the slope dM/dH and its scaling col- 
lapse. By using this derivative, the critical region is em- 
phasized as the peak in the curve, and the dependence 
on the parameter M c drops out. The lower graphs in 
figure ED and ^| show the scaling collapses of the mag- 
netization and its slope. Clearly in neither case is all the 
data collapsing onto a single curve. This would be dis- 
tressing, were it not for the fact that this also occurs in 



12 This section follows closely the presentation in J26 



FIG. 15 M(H) curves at different disorders, and (poor) 
universal scaling collapse. ( 28) 



mean field theory <|43r> at a similar distance to the critical 
point. 

Because the scaling of the magnetization is so bad, we 
use other quantities to estimate the critical exponents 
and the location of the critical point (Tables I and III) . 
Fixing these quantities, we use the collapse of the dM/dH 
curves to extract the rotation B mixing the experimental 
variables r and h into the scaling variable h! = h + Br 
( equation IA13I and following discussion). Recent work 
suggests that much of the difficulty in collapsing the mag- 
netization curves is surmountable by adding analytical 
corrections to scaling (|44T) . 



B. Avalanche Size Distribution 

In our model the spins flip in avalanches: each spin 
can kick over one or more neighbors in a cascade. 
These avalanches come in different sizes. The integrated 
avalanche size distribution is the size distribution of all 
the avalanches that occur in one branch of the hystere- 
sis loop (for H from — oo to oo). Figure [31 (3) shows 
some of the raw data (thick lines) in 3 dimensions. Note 
that the curves follow an approximate power law behav- 
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FIG. 16 ^-{H) curves at different disorders, and 
(poor) universal scaling collapse. l|2Si ) While the curves 
are not collapsing onto a single curve, the quality of the col- 
lapse is quite similar to that found at similar distances from 
R c in mean field theory J43fl . for which we know analytically 
that scaling works as R — ► R c . 



ior over several decades. Even 50% away from criticality 
(at R = 3.2), there are still two decades of scaling, which 
implies that the critical region is large. In experiments, 
a few decades of scaling could be interpreted in terms of 
self-organized criticality. However, our model and simu- 
lation suggest that several decades of power law scaling 
can still be present rather far from the critical point (note 
that the size of the critical region is non-universal) . The 
slope of the log- log avalanche size distribution at R c gives 
the critical exponent r + a/38. Notice, however, that the 
apparent slopes in figure [3] continue to change even after 
several decades of apparent scaling is obtained. The cut- 
off in the power law diverges as the critical disorder R c 
is approached. This cutoff size scales as S ~ |r|~ 1 / <T . 

These critical exponents can be obtained by using a 
scaling collapse for the curves of figure [31 shown in the 
inset. The scaling form for the avalanche sizes as a func- 
tion of R and H is 



We can find the scaling form for the distribution of 
avalanche sizes integrated over H by integrating this for- 
mula, changing variables to y = h/r" , and rewriting 
the resulting formula in terms of S to a power times a 
function of S a r. For R> R r , 



D in t(S,R) 



D{S, R, H) dH 



(10) 



oc / S- T V(S a r,h/r^)dh 



j-VM 



V(S"r,y)dy 






..(mi) 



where T> + is the scaling function for the integrated 
avalanche size distribution (the + sign indicates that the 
collapsed curves are for R > R c ). We are sufficiently 
far from the critical point that corrections to scaling are 
important: as described in reference l|43r) . we do collapses 
for small ranges of R and then linearly extrapolate the 
best-fit critical exponents to R c . We estimate from this 
curve that the critical exponents r + a/36 = 2.03 and 
a = 0.24 

The scaling function T>^ nt) (X) with X = S a \r\ is a 
universal prediction of our model. To facilitate compar- 
isons with experiments, we fit a curve to the data collapse 
in the inset of figure 01 We have fit the scaling collapses 
in dimensions 3, 4, and 5 to a phenomenological form of 
an exponential times a polynomial. In three dimensions, 
our fit is 



V ( ™ t] {X) = e 



-0.789X 1/o 



D(S, R, H) ~ S- T V{S a \rl h/r 



0S\ 



(9) 



(11) 
(0.021 + 0.002X + 0.531X 2 - 0.266X 3 + 0.261X 4 ) 

where 1/a = 4.20. The distribution curves obtained us- 
ing the above fit are plotted (thin lines in figureOU along- 
side the raw data (thick lines). They agree remarkably 
well even far above R c . We should recall though, that 
the fitted curve to the collapsed data can differ from the 
"real" scaling function even for large sizes and close to 
the critical disorder (in mean field Ij43j) the error in the 
corresponding curve was about 10%). The scaling func- 
tion in the inset of figure has a peculiar shape: it grows 
by a factor of ten before cutting off. The consequence of 
this bump in the shape is that in the simulations it takes 
many decades in the size distribution for the slope to 
converge to the asymptotic power law. This can be seen 
from the comparison between a straight line fit through 
the R = 2.25 (billion spin) simulation in figure |3 and the 
asymptotic power law S~ 203 obtained from extrapolat- 
ing the scaling collapses (thick dashed straight line in the 
same figure). A similar bump exists in other dimensions 
and mean field as well. FigureEl s h ows the scaling func- 
tions in different dimensions and in mean field. In this 
graph, the scaling functions are normalized to one and 
the peaks are aligned (the scaling forms allow this) . The 
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curves plotted in figure El are n °t raw data but fits to 
the scaling collapse in each dimensions, as was done in 
the inset of figure [3 
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FIG. 17 Universal scaling functions X'int(<S" T |r|) for the 
avalanche size distributions in different spatial dimen- 
sions. d28T) The inset of figure [3] shows the scaling collapse for 
three dimensions. 

It is clear from the figure that the growing bump in the 
scaling curves as the dimension decreases is a foreshad- 
owing of a zero in the scaling curve in two dimensions. 



C. Avalanche Correlations 
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The avalanche correlation function G(x, R, H) mea- 
sures the probability that the initial spin of an avalanche 
will trigger, in that avalanche, another spin a dis- 
tance x away. From the renormalization group descrip- 
tion l|27t J45f) . close to the critical point and for large 
distances x, the correlation function is given by: 



G(x,R,H) 



1 



d-2+i] 



X 



Q±{x/Z{r,h)) 



(12) 



where r and h are respectively the reduced disorder and 
field, Q± (± indicates the sign of r) is the scaling function, 
d is the dimension, £ is the correlation length, and r\ is 
called the "anomalous dimension" . Corrections can be 
shown to be subdominant ( appendix IX|l . The correlation 
length £(r, K) is a macroscopic length scale in the system 
which is on the order of the mean linear extent of the 
largest avalanches. At the critical field H c (h=0) and 
near R c , the correlation length scales like £ ~ \r\~" , while 
for small field h it is given by 



^ \rr y ± (h/\rf s ) 



(13) 



where y± is a universal scaling function. The avalanche 
correlation function should not be confused with the clus- 
ter or "spin-spin" correlation which measures the prob- 
ability that two spins a distance x away have the same 



FIG. 18 Avalanche correlation function integrated over 
field H at various disorders, with scaling collapse (2&). 



value. (The algebraic decay for this other, spin-spin cor- 
relation function at the critical point (r = and h = 0), 
is l/x d - 4+f > gl).) 

We've mostly used, for historical reasons, a slightly 
different avalanche correlation function, which scales the 
same way as the "triggered" correlation function G de- 
scribed above. Our function basically ignores the differ- 
ence between the triggering spin and the other spins in 
the avalanche: alternatively, it calculates for avalanches 
of size S the correlation function for pairs of spins, 
and then averages over all avalanches (weighting each 
avalanche equally). We've checked that our correlation 
function agrees to within 3% with the "triggered" cor- 
relation function described above, for R > R c in three 
dimensions and above. (In two dimensions, the two def- 
initions differ more substantially, but appear to scale in 
the same way.) 

We have measured the avalanche correlation function 
integrated over the field H, for R > R c . For every 
avalanche that occurs between H = — oo and H = +oo, 
we keep a count on the number of times a distance x 
occurs in the avalanche. The spanning avalanches are 
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not included in our correlation measurement. Figure [TBI 
shows several avalanche correlation curves in 3 dimen- 
sions for L — 320. The scaling form for the avalanche 
correlation function integrated over the field H, close to 
the critical point and for large distances x, is obtained 
by integrating equation Ijl2(l : 



G int (x,R) 



1 



:r 



d-2+r) 



G±(x/£(r,h)) dh (14) 



Using equation (|13|) and defining u 
(|T4*jl becomes: 



h/\r\" , equation 



G int (x,R) ~ \rf 5 x-( d - 2+ ^ I ' G±(x/\r\-»y±(u)) 



du 

(15) 

The integral (I) in equation (|15|) is a function of x\r\ v 
and can be written as: 



l=(x\rn 
to obtain the scaling form: 



-/38/v 



G±{M V ) 



G int (x,R) 



1 

x d+(3/v 



g±(x\r\n 



(16) 



(17) 



(35-13 



where we have used the scaling relation (2 — 77)1' 
(see (|45j) for the derivation). 

The bottom figure HH1 shows the integrated avalanche 
correlation curves collapse in 3 dimensions for L = 320 
and R > R c . The exponent v is obtained from such 
collapses by ex trap olating to R = R c as was done for 
other collapses (43). The exponent (3/v can be obtained 
from these collapses too, but we found it better to use 
the jump in magnetization l)28f) . which near the criti- 
cal point involves several spanning avalanches. Recent 
work pot 147c Bql suggests that the jump in magnetiza- 
tion may be dominated by only one of these spanning 
avalanches, and their work suggests that there may be 
two exponents related to our (3, one substantially larger 
(their (3 C ~ 0.15±0.08). The value of 0/v listed in Table I 
is derived exclusively from the magnetization discontinu- 
ity collapses. 

One of the most artificial features of our model is that 
the domains are arranged in a regular grid. It is thus 
important to check that this grid is not important for 
our scaling properties: if the avalanches looked notice- 
ably cubical, we would have problems. We have looked 
for possible anisotropies in the integrated avalanche cor- 
relation function in 2 and 3 dimensions. The anisotropic 
integrated avalanche correlation functions are measured 
along "generalized diagonals" : one along the three axis, 
the second along the six face diagonals, and the third 
along the four body diagonals. We compare the inte- 
grated avalanche correlation function and the anisotropic 
integrated avalanche correlation functions to each other, 
and find no anisotropies in the correlation, as can be seen 
from figure iTT?l Spherical symmetry emerges at the criti- 
cal point on long length and time scales. 




10 
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FIG. 19 Anisotropies in the avalanche correlation 
function. Four curves are shown on the graph: one is the 
avalanche correlation function integrated over the field H (as 
in figure 1181 , while the other three are measurements of the 
correlation along the three axis, the six face diagonals, and 
the four body diagonals. Avalanches involving more than four 
spins show no noticeable anisotropy: the critical point ap- 
pears to have spherical symmetry. The same result is found 
in 2 dimensions. 



D. Avalanche Duration Measurement 

The scaling relation between the duration T of an 
avalanche (the time it takes to occur) and the linear size 
£ of the avalanche defines the critical exponent z 



T~£ 2 



(18) 



The exponent z is known as the dynamical critical ex- 
ponent. Equation i|18fl gives the scaling for the time it 
takes for a spin to "feel" the effect of another a distance £ 
away. Since the correlation length £ scales like r~ v close 
to the critical disorder, and the characteristic size S as 
r -1 /' 7 , the duration T then scales with avalanche size as: 



T~S° 



(19) 



In our simulation, we measure the distribution of 
durations for each avalanche size S. The probability 
D t (S, R, H, T) that an avalanche of size S will be of du- 
ration T close to the critical field H c and critical disorder 
R c has the scaling form 



D t {S, R, H, T) ~ S- q V ( l\s°\rl h/\rf s , T/S° 
where q = r + uvz, and is defined such that 

/ + OO />OC 
/ D t (S,R,H,T) dH dt = 
-00 J\ 



(20) 



s -(T+*f38) D^*)(^| r |) 



(21) 



■*(int) 



where T>± was defined in the integrated avalanche size 
distribution section. The avalanche time distribution in- 
tegrated over the field H, at the critical disorder (r = 0) 
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D 



(int) 



/c r r\ ro j.— (T+cr/38+cruz)/auz j\\ m V Irnj qovz\ (00\ 



as obtained from equation 12U|) (reference l|43|) '). 




FIG. 20 Distribution of avalanche durations for 
avalanches of fixed size 5, for various sizes 5*, and 
their scaling collapse. (28) 

Figures 1201 show the avalanche time distribution inte- 
grated over the field H for different avalanche sizes, and 
a collapse of these curves using the above scaling form. 



E. Energy Spectrum 

The energy spectrum of the voltage V(t) as a function 
of time for our model (given by adding up the the voltage 
traces of each avalanche, as in figure |SJ), 



E(u) 



■(V(t)V(t + T)) = \V(u)\ 



(23) 



is a commonly measured experimental quantity 13 There 
are two distinct contributions to this power spectrum: 




FIG. 21 Power spectrum scaling and relatives. Top: 
short-range model; bottom: mean field (long-range dipolar 
forces (113).) 



(1) at high frequencies, the incoherent sum of the power 
spectra E mc of the individual avalanches, and (2) at low 
frequencies, a term representing the correlations between 
avalanches. As the system is forced more and more 
slowly, the internal dynamics of the avalanches is un- 
changed but their separation in time increases, so these 
two contributions separate. Typical experiments seem to 
be dominated by the contribution (1) from the internal 
dynamics within an avalanche (|13J) . In this section, the 
results we discuss are valid for all three models of hystere- 
sis. 14 See also (42) in this volume for a more complete 
discussion of our power spectrum theory and the experi- 
mental context. 



13 The power spectrum is E(uj) divided by the integration time. 



14 In the short— range model, we concentrate on the case where we 
measure only near R c and H c : we do not integrate over the loop. 
Integrating over the loop yields r very close to two, where both 
forms of the scaling should compete. 
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There have been a number of naive arguments for the 

ower law for E(ui), stretching back to 1972 H IH III 

49), which yield a power law which is only valid if r > 2: 



E u 



,(w) oc w -(3-t)/«'* ( only for T > 2 ). (24) 



Let ^max be the typical largest avalanche in the system, 
say because the system is finite in size. At root, these 
derivations fail for t < 2 because 



s„ 



SV(S) dS = 



S 1 ^ dS = Si 



(25) 



diverges. Since there are S spins in each avalanche, it 
would seem that equation|2Slshould represent the number 
of spins in the system. This is roughly correct for r > 2, 
but for t < 2 the amplitude of the power law D(S) = 
DqS~ t depends on this cutoff Dq ~ 5 m ^ x . Continuing a 
detailed analysis l|13r . we get a different scaling form 



E„ 



,(oj) <x 0J- 1/avz (forr<2). 



(26) 



It is not a surprise that this new prediction fixes a number 
of discrepancies between theory and experiments fl3|) . 

Figure |^ shows the power spectrum P(u>) (equal 
to E(oS) divided by the duration of the measurement), 
together with four other curves, for both the front 
propagation model and the infinite-range model l|13r). 
The top two curves are the naive power law (eauation l24l 
obviously not correct) and the correct law ( equation I26fl . 
We compare with the relation between the avalanche size 
S and the duration T feauation ll9[) . which has the same 
combination of critical exponents S ~ (l/T) _1 / CT,yz . We 
also show that the energy spectrum of individual large 
avalanches is proportional to S and has the same power 
law in frequency to as that of the entire time series, so we 
plot P(u)\S)/S for a variety of avalanche sizes S. 



F. Tables of Results 



Here we summarize, from l|28fl . our numerical estimates of the universal critical exponents and various non-universal 
quantities used in the scaling collapses above. 15 

measured exponents 3d 4d 5d mean field 



\/u 


0.71 ±0.09 


1.12±0.11 


1.47±0.15 


2 




6 


0.015±0.015 


0.32±0.06 


1.03 ±0.10 


1 




(t + a/38 - 3)/av 


-2.90 ±0.16 


-3.20 ±0.24 


-2.95 ±0.13 


-3 




1/a 


4.2±0.3 


3.20±0.25 


2.35±0.25 


2 




t + <t(38 


2.03±0.03 


2.07±0.03 


2.15±0.04 


9/4 




T 


1.60 ±0.06 


1.53 ±0.08 


1.48 ±0.10 


3/2 




d + fi/v 


3.07±0.30 


4.15±0.20 


5.1 ±0.4 


7 (at d c = 


= 6) 


p/v 


0.025 ±0.020 


0.19±0.05 


0.37 ±0.08 


1 




avz 


0.57 ±0.03 


0.56 ±0.03 0.545 ±0.025 


1/2 





Table I. Measured universal critical exponents. Values for the exponents extracted from scaling collapses in 3, 4, 
and 5 dimensions. The mean field values are calculated analytically (23; 45). v is the correlation length exponent and 
is found from collapses of avalanche correlations, number of spanning avalanches, and moments of the avalanche size 
distribution data. The exponent 8 is a measure of the number of spanning avalanches and is obtained from collapses 
of that data, (r + a (38 — i)/av is obtained from the second moments of the avalanche size distribution collapses. 
1/(7 is associated with the cutoff in the power law distribution of avalanche sizes integrated over the field H , while 
t + (7/36 gives the slope of that distribution, r is obtained from the binned avalanche size distribution collapses (|28T) . 
d + P/v is obtained from avalanche correlation collapses and 0/v from magnetization discontinuity collapses, avz is 
the exponent combination for the time distribution of avalanche sizes and is extracted from that data. Error bars are 
based on variations in the results based on different approaches to the analysis: statistical fluctuations are typically 
smaller. 



15 Some of the scaling collapses shown use d exponents slightly dif- 
ferent from those in the tables: see l28t). 
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calculated exponents 



3d 



4d 



5d 



mean field 



cr/3(5 
05 


(TV 



0.43 ±0.07 0.54 ±0.08 0.67 ±0.11 

1.81 ±0.32 1.73±0.29 1.57±0.31 

0.035 ±0.028 0.169 ±0.048 0.252 ±0.060 

0.34 ±0.05 0.28 ±0.04 0.29 ±0.04 

r) = 2 + (0-06)/p 0.73±0.28 0.25±0.38 0.06±0.51 
Table II. Calculated universal critical exponents. Values for exponents in 3, 4, and 5 dimensions that are 
not extracted directly from scaling collapses, but instead are derived from Table I and the exponent relations (see 
section IVH . The mean field values are obtained analytically f23t l45|) . Both a 05 and (3d could have larger systematic 
errors than the errors listed here <|28|) . Recent work l[4a) on small systems, but with sophisticated scaling analysis, 
suggests a larger value for (3. 



3/4 
3/2 
1/2 
1/4 




3d 



4d 



5d 



field 



R c 2.16 ±0.03 4.10 ±0.02 5.96 ±0.02 0.79788456 

H c 1.435 ±0.004 1.265 ±0.007 1.175 ±0.004 

B 0.39 ±0.08 0.46 ±0.05 0.23 ±0.08 

Table III. Non-universal scaling variables. Numerical values for the critical disorders and fields, and the rotation 
parameter B ( equation IA13I and subsequent discussion), in 3, 4, and 5 dimensions extracted from scaling collapses. 
The critical disorder is obtained from collapses of the spanning avalanches and the second moments of the avalanche 
size distribution. The critical field is obtained from the binned avalanche size distribution ( 28) and the magnetization 
curves. H c is affected by finite sizes, and systematic errors could be larger than the ones listed here. The mean field 
values are calculated analytically <|23tl45l) . The rotation B is obtained from the dM/dH collapses. 



r 



VI. EXPONENT RELATIONS 

In the following sections we list various exponent re- 
lations for the nonequilibrium RFIM, for which we give 
detailed arguments in <|27t 1501) . 



A. Exponent equalities 

The exponents introduced above are related by the fol- 
lowing exponent equalities: 



0-/36= (r-2)/crif r< 2, 



(2 - rj)v = 05 -(3, 



(27) 



(28) 



B. Two violations of hyperscaling 

In the nonequilibrium RFIM there are two different 
violations of hyperscaling. 

(1) In l)27t I50T). we show that the connectivity hyper- 
scaling relation 1/er = dv — from percolation is violated 
in our system. There is a new exponent 9 defined by 



1/(7= (d-0)v- 



(31) 



with Ov = 1/2 - e/6 + 0(e 2 ) and 9v = 0.021 ± 0.021 in 
three dimensions. 9 is related to the number of system 
spanning avalanches observed during a sweep through the 
hysteresis loop: see also Ma). 

(2) As we discuss in (|27l ) there is a mapping of the 
perturbation theory for our problem to that of the equi- 
librium random-field Ising model to all orders in e. From 
that mapping we deduce the breakdown of an infamous 
( "energy" )-hyperscaling relation 



P=\^ 



n), 



and 



(29) 



5 = (d - 1r\ + fj)/{d -i + fj). 



(30) 



(The latter three equations are not independent and 
are also valid in the equilibrium random-field Ising 
model) . 



+ 05=(d- 0)u, 



(32) 



with a new exponent 9, which has caused much contro- 
versy in the case of the equilibrium random-field Ising 
model. 

In Ij27t 15a) we discuss the relation of the exponent 9 
to the energy output of the avalanches. The e expansion 
yields 9 = 2 to all orders in e. Nonperturbative correc- 
tions are expected to lead to deviations of 9 from 2 as 
the dimension is lowered. The same is true in the case 
of the equilibrium RFIM. The numerical result in three 
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dimensions is 9 = 1.69 ± 0.28. (In the three-dimensional 



equilibrium RFIM ||5J) it is eq = 1.5 ± 0.4). 

Another strictly perturbative exponent equality, which 
is also obtained from the perturbative mapping to the 
random-field Ising model is given by 



77 = 77. 



(33) 



It, too, is expected to be violated by nonperturbative 
corrections below six dimensions. 



C. Exponent inequalities 



In Ij27t |50J) we give arguments for the following two 
exponent-inequalities. 16 



v/(35 > 2/d, 



(34) 



which is formally equivalent to the "Schwartz-Soffer" 
inequality, v < 2u, first derived for the equilibrium 
random-field Ising model, and 



the RFIM), then at (small enough) finite field sweep rate 
the corresponding noise "pulse" duration distribution is 
expected to scale as D(T, fi) ~ T 2 - a ^ cn where a(H) 
and C are nonuniversal constants l|54t l55t l56t l57j) . 

If, however, in the adiabatic limit the exponent a is 
either greater or smaller than 2, then at (small) finite 
field sweep rate the pulse duration distribution exponent 
remains the same as in the adiabatic limit. Note that in 
the case a < 2 especially, f2 has to be particularly small in 
order to be able to still see distinct pulses even at finite 
field sweep rate, since at higher sweep rate the system 
will very quickly develop runaway events. In 11521 : |53J) 
quantitative criteria for the meaning of "small sweeprate" 
are given. Also, the zero temperature nonequilibrium 
RFIM and recent variants are used to numerically test 
the analytic results, which are expected to be applicable 
to a much larger class of systems with crackling noise 
(the exact conditions are discussed in Ij52t l53m. It is 
found that the results agree well with both simulations 
and recent experiments on Barkhausen noise in various 
soft magnetic materials with and without applied stress 
(with a < 2 and a — 2 respectively iflllO. A brief review 
of other, related studies of the effects of finite field sweep 
rate is also given in l|52r) . 



v > 2/d, 



(35) 



which is a weaker bound than Eq. (3) so long as (35 > 1, 
as appears to be the case both theoretically and numer- 
ically at least for d > 3. 



VII. FINITE SWEEPRATE 



VIM. SUBLOOPS AND HISTORY INDUCED CRITICAL 
BEHAVIOR 

One of the characteristic features of magnetic hystere- 
sis are the subloops, seen as the external field is changed 
up and down an amount insufficient to saturate the mag- 
netization (figure 0J. 



Originally the nonequilibrium RFIM was studied in the 
adiabatic limit, where the external field is kept constant 
during each spin flip avalanche and only increased be- 
tween subsequent avalanches. In real experiments, how- 
ever, the driving field H is typically increased at a finite 
rate fi such that H = fit where t is time. This finite 
rate allows for new avalanches to be triggered by the in- 
creasing external field while an earlier avalanche is still 
propagating. In Ij52t |53J) the effect of finite field sweep 
rates on the power spectra and avalanche size distribu- 
tions is discussed in detail for a large class of systems with 
crackling noise. In particular, it is asked how the scaling 
behavior of the avalanche size and duration distribution 
and the power spectra of Barkhausen noise or crackling 
noise in general depend on the field sweep rate fl. One of 
the results is an exponent inequality as a criteria for the 
relevance of adding a small driving rate fl > to the adi- 
abatic case O — > 0: If in the adiabatic case the avalanche 
duration distribution scales as D(T) ~ T~ Q with a = 2 
(such as in the ABBM model or the mean field version of 



16 From the normalization of the avalanche size distribution 
D(s, r, h) it follows that r > 1. 



A. Return-Point Memory 

Much attention has been paid to modeling these 
subloops with the Preisach model IJ58T) Preisach models 
are quite different in spirit to ours: they represent the 
system as a large number of uncoupled hysteretic two- 
state domains, and fit the distribution to the observed 
behavior. They are able to model a hysteretic system 
if it possesses return-point memory IjlaK also known as 
wiping out l|59|) . 

Return-point memory states that the magnetization 
after a subloop rejoins a larger loop (perhaps itself a 
subloop) equals the magnetization at which the subloop 
left the outer loop. That is, if the subloop represents 
an excursion downward from a local maximum exter- 
nal field H, then when H returns to its previous max- 
imum the magnetization M returns to its previous value 
- remembering its previous magnetization on returning, 
and wiping out all effects of the excursion. Return-point 
memory states that the subloops in figure 01 should close 
perfectly, without a gap and without crossing themselves 
(as observed). 

Indeed, our interacting, three-dimensional, disordered 
models exhibit the return-point memory in an even 
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stronger form. Upon rejoining the larger loop, the en- 
tire state of the system is microscopically identical to the 
state it had when it left the larger loop (a microscopic 
return-point memory, as opposed to a macroscopic mem- 
ory constraining only the net magnetization). In l|23f) . 
we showed in great generality that return point memory 
should be true of any system with 

1. A partial ordering of states. Here, one mi- 

crostate A is "ahead" of another B, A > B, if every 
spin in A is greater than or equal to the correspond- 
ing spin in B. 

2. No passing. That is, the dynamics preserves the 
partial ordering. 

3. Adiabatic. The external field is raised and low- 
ered slowly enough that the system does not lag 
behind. 



—Hr, in the form of a kink in the reversal of the mag- 
netization of the reversal curve. They show that this 
memory effect is due to a local spin-reversal symmetry of 
the Hamiltonian. This symmetry and thus the reversal 
field memory effect is present in spin glasses where the 
disorder is due to random couplings between spins. In 
general it is not expected in systems with random mag- 
netic fields such as the RFIM, since there the disorder 
breaks local spin reversal symmetry. 

As an extension to memory effects for a driving field 
that varies in two dimensions, a vector form of wipe-out 
memory for situations when the magnetic field varies in 
both direction and magnitude, has been suggested based 
on the two dimensional vector Preisach model in (72 
and has been experimentally tested in 1)731) . 



B. Critical Behavior in Subloops 



Return-point memory can be quite remarkable: in our 
model (and in some experimental systems) repeating a 
subloop plays back precisely the same Barkhausen noise. 
Indeed, the return-point memory can be used as a way of 
retrieving analogue magnetic memories that is (at least 
theoretically) significantly superior to measuring the re- 
manent magnetization (60). On the other hand, the ab- 
sence of exact (microscopic) return point memory has 
also been observed (|61t l62t l63j) . and there are various 
interesting experiments testing for the reproducibility of 
magnetic avalanches IjrM l65i l66f) and for microscopic re- 
turn point memory and complementary point memory 
in magnets at various disorders l)67|) . 17 In (J67|) the au- 
thors experimentally study the influence of disorder on 
major loop return point memory and complementary re- 
turn point memory in Co/Pt samples with varying in- 
terfacial roughness and find with increasing disorder the 
onset and saturation of both return point memory and 
complementary point memory. 

Disorder dependence of return point memory versus 
reptation (i.e. is gradual subloop closure (68; 69)) has 
recently been reported for the zero temperature ran- 
dom coercivity model with antiferromagnetic-like inter- 
actions 1|70|) . 

A related memory effect, recently reported for the hys- 
teresis of spin glasses, is the reversal field memory ef- 
fect (I71f) . Katzgraber et al. show that this memory effect 
emerges in the nonequilibrium Edwards Anderson spin 
glass (EASG) when the magnetic field is first decreased 
from its saturation value and then increased again from 
some reversal field Hr. The authors find that EASG ex- 
hibits a singularity at the negative of the reversal field, 




-0.5 



FIG. 22 Saturation hysteresis loops with concentric 
inner subloops. Top: R = 2.225J > R c ; bottom: R = 
2.1 < R c . (74) 



17 The complementary point memory relates the magnetic domains 
at one point on the major hysteresis loop to the domains at the 
complementary point on the major loop during the same and 
during subsequent cycles. 



The disorder induced critical scaling for hysteresis 
loops in the RFIM that we discussed so far all referred 
to the saturation hysteresis loop. However, often it is 
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experimentally impractical to take magnets to their sat- 
uration point due to the large magnetic fields required, 
so the behavior of subloops (figure E3 |Z3)) is 0I great 
interest to experiments and applications. The RFIM can 
be used to model subloops, and in one dimension l|7q) 
and on a Bcthc lattice l|7rl Yur i they have been com- 
puted exactly l|78|) . The magnetization curves of subloops 
have been collapsed near the demagnetized state using 
Rayleigh's law {SHU)- Interestingly, in the RFIM not 
only the saturation hysteresis loops but also subloops re- 
flect the critical point described above, and at the crit- 
ical disorder one finds history-induced critical scaling. 
Concentric inner subloops resemble rescaled saturation 
loops at effectively higher (possibly correlated) disorder. 
Avalanche size distributions for the inner subloops, too, 
look similar to rescaled avalanche size distributions of 
saturation loops at higher disorder. We have studied 
universal properties of this history induced critical be- 
havior using numerical simulations of systems with more 
than 100,000,000 spins and Widom scaling methods and 
extracted universal scaling exponents and scaling func- 
tions. These predictions are expected to apply to ex- 
periments on hysteresis and Barkhausen noise in hard 
magnetic materials. To model soft magnets on the other 
hand, we added demagnetizing fields to the model and 
obtained single domain wall propagation, as described 
in (13). As a result the system was always critical, even 
for the subloops. The results of this study are discussed 
in <|74l: l80|) . A special, related kind of subloops - so called 
"first order reversal curves" (FORC) - have proven very 
useful in distinguishing between different magnetic ma- 
terials, see Jill). 



C. Demagnetization Curves 

The RFIM also predicts disorder induced critical be- 
havior in the demagnetization curve, or initial magneti- 
zation curve and the associated Barkhausen noise. Fer- 
romagnetic materials with a remanent magnetization can 
be demagnetized by applying an oscillating magnetic field 
with amplitude slowly decreasing from a large initial 
value to zero. Sometimes the final state is also termed 
AC-demagnetized state. The oscillating external field 
with decreasing amplitude takes the system through con- 
centric subloops. The line connecting the tips of the 
subloops is known as the normal or initial magnetiza- 
tion curve, or demagnetization curve. It is to be dis- 
tinguished from the virgin curve which is obtained by 
thermal demagnetization q 16c l82t |83j). One can use the 
RFIM to simulate the demagnetization curve for differ- 
ent amounts of disorder, especially near the critical dis- 
order mentioned above for the saturation loop. It turns 
out that this curve reflects much of the scaling behavior 
found for the saturation hysteresis loop. This may be 
surprising, since the meta-stable states encountered in 
the demagnetization curve are completely different from 
those of the saturation loop. In (841 ) the model is briefly 



described and the necessary exponent relations are de- 
rived, which show that the exponents for the demagne- 
tization curve can actually be derived from those of the 
subloops. Predictions for the exponents from simulations 
of up to one million spins arc presented. The demagne- 
tization curve in the presence of long range demagnetiz- 
ing fields is also discussed, with resulting predictions for 
experiments on demagnetization curves and associated 
Barkhausen noise in soft magnetic materials. For a dis- 
cussions of the Rayleigh law and scaling of the demag- 
netization curve at small magnetic fields, see l|77t J79f) . 
Zarand, Pazmandi, Pal, and Zimanyi recently proposed 
a new optimization method to search for ground states, 
which is based on the demagnetization procedure (|85r) - 



IX. REAL EXPERIMENTS 
A. Magnetic hysteresis loops 

There are a few recent experiments on hysteresis and 
crackling noise where disorder is tuned and evidence of 
a phase transition related to that of the nonequilibrium 
RFIM is reported: In (86), Berger et al. describe a tran- 
sition from smooth magnetic hysteresis loops to those 
with a jump in the magnetization in ultrathin epitaxial 
Gd-films with in-plane magnetization upon annealing the 
sample at successively higher annealing temperatures. 
The authors show that annealing at higher temperature 
increasingly reduces the disorder in the crystal structure 
of the films. In lJ87t l88t l89f) Berger et al. report a simi- 
lar transition for thin Co/CoO films. In this system the 
antiferromagnetic CoO layer allows a reversible tuning of 
the magnetic disorder felt by the ferromagnetic Co layer, 
by simple temperature variation. The authors are able 
to extract scaling collapses for the magnetization hys- 
teresis loops at different disorders, and find the values of 
the associated scaling exponents to be f3 — 0.022 ± 0.006 
and (35 = 0.30 ± 0.03. Note that this magnetic system is 
effectively two dimensional. 

Incidentally, almost the same exponent values (j3 — 
0.03 ± 0.01 and fid = 0.4 ± 0.1) are reported by Marcos 
et al. in I|9CJ1 for the magnetic hysteresis loops of Cu-Al- 
Mn alloys with different Mn content at low temperatures. 
The loops are smooth above a certain temperature, but 
exhibit a jump in the magnetization below that temper- 
ature. The authors perform a scaling analysis near this 
critical point with temperature as the effective "disorder" 
tuning parameter. They point out that, contrary to the 
thin Co/CoO bilayers in Berger et al. <|87J) . this system 
is effectively three dimensional. 

For these three systems it would be very interesting to 
also see a scaling analysis of the associated Barkhausen 
noise, since from simulations we know that the expo- 
nents extracted from avalanche size distributions and 
noise power spectra are often more reliable than those 
obtained just from magnetization curves alone IJ27J) . For 
these three systems, so far, only the scaling of the mag- 
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netization curves has been reported. 



for thin films with in-plane applied magnetic field, or for 
thin wires with parallel applied field). 



B. Disorder effects on Barkhausen Noise 

Disorder effects on crackling noise (rather than hys- 
teresis loops) has been seen in other systems. Barkhausen 
noise is commonly used for nondestructive testing - in 
engineering contexts we find careful studies of the effect 
of the amount of carbon content in steel on Barkhausen 
noise. One might speculate that varying the carbon con- 
tent could be similar to varying the effective disorder in 
the system. It would be interesting to see a systematic 
scaling analysis for the available experimental data. In an 
early study IJ49T) Lieneweg and Grosse-Nobis found that 
the distribution of Barkhausen pulse areas integrated 
over the hysteresis loop of an 81% Ni-Fe wire was well 
described by a power law (with exponent r + aj35 ranging 
from 1.73 to 2.1) up to a certain cutoff size. Annealing 
the sample at various annealing temperatures they found 
that the cutoff appeared to be smaller at higher anneal- 
ing temperatures. It would be interesting to see whether 
the cutoff grows to takes a (system-size dependent) maxi- 
mum value at a critical annealing temperature T^ nn and 
decreases again at higher and lower annealing temper- 
atures, as one would expect if the RFIM applies with 
temperature as the effective disorder tuning parameter in 
this system. (In leading approximation, near T" nn they 
should be linearly related to each other.) Near T° n ™ the 
Barkhausen pulse area distribution should then be de- 
scribed by a scaling form that would allow a scaling col- 
lapse of all distributions onto one single curve for appro- 
priate rescaling of the axes. Again, potentially universal 
critical exponents could be extracted from such a col- 
lapse. They would be predicted by the RFIM avalanche 
critical exponents if that model is in the same universal- 
ity class. 

There are other experiments which revealed power law 
decays for Barkhausen pulse size distributions in various 
samples, as reviewed for example in OJ) and in the article 
by Durin and Zapperi in this book (J42J). Many of these 
are measurements in soft magnetic materials where long 
range interactions are important and the main mecha- 
nism for Barkhausen noise is domain wall propagation 
rather than new domain nucleation. Disorder has usu- 
ally not been varied in these experiments. Models pre- 
dict that these systems naturally tend to operate close to 
the underlying critical domain wall depinning transition 
and thus produce power law pulse size distributions re- 
gardless of the underlying disorder, (at least as long as it 
is less than a critical value above which the domain wall 
may have overhangs H 5 IS 113 E2 E|) The disor- 
der induced critical point of the noncquilibrium RFIM is 
expected to be relevant for hard magnetic materials in- 
stead, where either the disorder is very strong compared 
to the effects of long range dipolar interactions, or the 
geometry of the sample and experiment is chosen such 
that the effects of long range dipolar fields are small (as 



C. Imaging magnetic avalanches and states 

The avalanche structures predicted by the theoretical 
models are visually interesting f figure |2*3*|) : in our model 
they have fractal dimensions 1/av close to, but proba- 
bly less than three. Exciting experiments probing these 
spatial structures are coming on line. 




FIG. 23 Fractal spatial structure of an avalanche. This 
moderate-sized avalanche contains 282,785 domains (0). 



Recently interesting experiments imaging avalanches 
in magnetic films have been reported. In (91) Kim, 
Choe, and Shin report full-field magneto-optical obser- 
vations (via Kerr effect) of Barkhausen avalanches in 
Co polycrystalline thin films at criticality, i.e. their im- 
aged avalanches are due to single domain wall propaga- 
tion rather than domain nucleation. In this experiment 
many avalanches extend far beyond the imaged region, so 
that their measured exponent appears to be an exponent 
combination of the exponents computed in the models. 
In addition to their results it'd be interesting to also see 
time resolved images of entire avalanches and direct mea- 
surements of the roughness exponent of the propagating 
domain wall, to compare with predictions from models 
for magnetic domain wall propagation. The magneti- 
zation reversal dynamics of epitaxial Fe films grown on 
GaAs(OOl) and the field dependent evolution of domain 
structure at various field has also been studied using the 
magneto-optic Kerr effect, see for example l|92l) . In (|93|) 
Schwarz et al. report on visualization of the Barkhausen 
effect by magnetic force microscopy for a granular thin 
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film with perpendicular anisotropy. The authors are able 
to localize individual Barkhausen volumes and distin- 
guish them as either newly nucleated or grown by domain 
wall propagation. They find a Gaussian size distribution 
of nucleated Barkhausen volumes indicating an uncorre- 
cted random process, and a power law distribution for 
grown Barkhausen volumes, reflecting critical behavior 
during domain wall motion. Another very interes ting ex- 
perimental approach is x-ray speckle metrology IJ67T) . In 
a recent preprint Pierce et al use this method to study 
correlations between microstates with different histories 
and test for microscopic return point memory of subloops 
in magnetic systems. 



D. Random Bonds and Random Anisotropy 

Renormalization group treatments, symmetry argu- 
ments (I27l) . and simulations I|94il95r) show, that the hys- 
teretic, zero temperature random bond Ising model with 
ferromagnetic mean bond strength is expected to be in 
the same universality class as the nonequilibrium RFIM. 
Vives et al. also point out that their numerical scaling 
exponents in two dimensions of the two models also have 
overlapping error bars with the corresponding exponents 
of Blumc-Emcry Griffiths models with disorder, in which 
the spins take three values ±1, (95). An analysis of the 
short time dynamics of first-order phase transitions in a 
disordered system, using the random bond Ising model 
is given in (|9q ) . The diluted ferromagnetic random Ising 
model (]97u has both nonuniform bond strengths that re- 
flect dilution due to nonmagnetic impurities, and random 
magnetic fields. Symmetry arguments suggest (|27l) that a 
model with uncorrelated random (ferromagnetic) bonds 
and random fields should again be in the same univer- 
sality class as the RFIM, since it has no new symmetries 
compared with the RFIM. Similarly, daSilveira and Kar- 
dar developed a O(N) vector spin model with random 
fields and showed through a renormalization group cal- 
culation that, in general, it can displays the same con- 
tinuous phase transition in the hysteresis loop as the 
RFIM f98), i.e. the two critical points are in the same 
universality class. 18 The reason is that the applied ex- 
ternal field breaks any added rotational symmetry l|27f) 
so that the behavior on long length scales is determined 
by the same the same symmetries in both models. Sim- 
ilarly the random anisotropy model of <|99j) . was shown 
numerically to be in the same universality class as the 
RFIM: this is again as expected, since the external ap- 
plied field not only breaks rotational symmetry but also 
time-reversal symmetry. Perhaps somewhat surprisingly, 
a transition from hysteresis loops with a jump to those 
without, like the one of the RFIM, is also seen inaplas- 
tic depinning model for charge density waves (|XO0T> . If 



They also find a new universality class, see section IXI 



in the random field Ising model a certain fraction of the 
bonds between neighboring spins is enhanced to a very 
large value, one obtains a simple model for exchange bias 

iHol . 

Hysteresis loops observed for the Ising glass with zero 
ferroma gnetic me an bond strength are smooth and con- 
tinuous l)l02l Il03[) - qualitatively this agrees with exper- 
imenta l hy steresis loops seen for example in AuFe sys- 
tems l)l04|) . For Ising glasses with strong ferromagnetic 
mean bond strength (compared to the width of the distri- 
bution of bond strengths, which is the disorder), a jump 
in the magnetization curve appears, just as in the low dis- 
order phase of the RFIM. The hysteresis loops resemble 
those seen for more concentrated AuFe syst ems ()104|) as 
well as CuMn and AgMn systems l|l05Ul06J) . A brief re- 
view of these and other glassy spin systems such as Ising 
spin glasses with long-range interactions and Heisenberg 
spin glasses with long-range interaction s in both theory 
and experiment is given for example in (103). 

An interesting hysteresis model with rather different 
properties than ours based on an infinite-range spin glass 
(Sherrington-Kirkpatrick model) is reported to show self- 
organized criticality; see l|l02j) for details. 



E. Other dynamics at high sweep rates 

Much interesting work has been done in both the- 
ory and experiment on dynamic hysteresis loops at fi- 
nite field sweep rate in the presence of temperature. At 
field sweep rates higher than a critical value that de- 
pends on sweep amplitude and temperature, the system 
transitions to a symmetry broken phase with a flatter 
hysteresis loop, centered around a nonzero average mag- 
netization. Similar dynamic transitions also are seen 
for pulsed and stochastically varying fields. Chakrabarti 
and Acharyya published a review arti cle on both models 
and experiments in this regime (|107T ). Some of the ad- 
dressed questions deal with scaling of the hysteresis loop 
area with field sweep frequency, amplitude, and temper- 
ature. While in the initial approaches to these issues 
quenched randomness in the system was neglected, there 
are recent model approaches and experiments studying 
dynamic hysteresis of a moving domain wall in the pres- 
ence of quenched disorder and temperature at finite field 
sweep rate ljl9il2l|) . 



F. Nonmagnetic noisy hysteretic systems 

In this section we briefly discuss a broad range of re- 
lated work on avalanches and hysteresis in non-magnetic 
systems. We start with a general overview, and then 
focus on a few particular systems. 
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1. Systems with avalanches 

Magnets crackle (Barkhausen noise) when the exter- 
nal field is ramped upward slowly. Many systems exhibit 
similar crackling noise: when pushed slowly they respond 
with discr ete e vents of a broad range of sizes. 19 The earth 
responds ( 108) with violent and intermittent earthquakes 
as two tectonic plates rub past one another. A p i ece o f 
paper Q) (or a candy wrapper at the movies (|l09t IllCf) ^ 
emits intermittent, sharp noises as it is slowly crumpled 
or rumpled. Just like the magnetization avalanches re- 
sponsible for Barkhausen noise in magnets, these indi- 
vidual events span many orders of magnitude in size. In- 
deed, the distribution of sizes follows a power law with 
no characteristic size scale. 

In the past decades, scientists in a variety of fields 
have been making rapid progress in developing models 
and theories for understanding this sort of scale- invariant 
behavior in driven, nonlinear, dynamical systems. In- 
terest in these sorts of phenomena goes b ack several 
decades. The work of Gutenberg and Richter IJ108I ) in the 
1940's and 1950's established the well-known frequency- 
magnitude relationship for earthquakes that bears their 
names. A variety of many-degree-of - freedom dynamical 

mod els iTTrTT: ITTa ITT1 [Til [Til ITT1 Imt ITT1 lul IT20I 

Il2ll Il22|) . with and without disorder, have been intro- 
duced in the years since to investigate the nature of slip 
complexity in earthquakes. More recent impetus for work 
in this field came from the study of the depinning tran- 
sition in sliding charge-density wave (CDW) c onductors 
in the 1980's and early 1990's (EH EH EH EH EH 
Il28t Il29j) . Interpretation of the CDW depinning tran- 
sition as a dynamic crit i cal p henomenon sprung from 
Fisher's early work l)l23t Il24|) . and several theoretical 
and numerical studies followed. This activity culminated 
in the RG solution by Narayan a nd Fi sher l|126r) an d the 
numerical studies by Middleton l)l27t) and Myers |)128|) 
which combined to provide a clear picture of depinning 
in CDWs and open the doors to the study of other dis- 
ordered, nonequilibrium systems. 

Bak, Tang, and Wiesenfeld ins pired much of the suc- 
ceeding work on crackling noise l|130i Il31|) . They in- 
troduced the connection between dynamical critical phe- 
nomena and crackling noise, and they emphasized how 
systems may naturally end up at the critical point 
through a process of self-organized criticality. (Their 
original model was that of avalanches in growing sand- 
piles. Sand has long been used as an example of crackling- 
noise l|l32Hl33j) . However, it turns out that r eal sandpiles 
don't crackle at the longest scales i|l34Ul35|) '). 

Researchers have studied many systems that crackle. 
Simple models have been developed to study bubbles re- 
arranging in foams as they are sheared IJ136D . biological 



extinctions l)l37|) , 20 fluids invading porous materials and 
other pro blems involving invading fronts Q l34t 1351 : E2 
141: ll42p . 21 the dyna mics of sup erconductors l|l43l Il44i 
145 ). and superfluids l)l46t Il47[). sound emitted during 
martensitic p hase transition s l|14cf) . fluct uatio ns in the 
stock market (|l49t Il50t Il5ll). solar flare s l|l52t ). cascad- 
ing failures in power grids !J153Hl54 : ll5.'il). failures i n sys - 
tems designed for optim al pe rformance l)l56t Il57t Il58f) . 
group decisio n ma king l)l59T) . crackling noise in mam- 
mali an lungs 11 1601) , and frac t ure i n disordered materi- 
als HH2 Il6l Il6fl EH3 EH EEB)- These models are 
driven systems with many degrees of freedom, which re- 
spond to the driving in a series of discrete avalanches 
spanning a broad range of scales what we are calling 
crackling noise. 

There has been healthy skepticism by some established 
professionals in these fields to the sometimes grandiose 
claims by newcomers proselytizing for an overarching 
paradigm. But often confusion arises because of the 
unusual kind of predictions the new methods provide. 
If our models apply at all to a physical system, they 
should be able to predict all behavior on long length and 
time scales, independent of many microscopic details f 
the real world. This predictive capacity comes, however, 
at a price: our models typically don't make clear predic- 
tions of how the real-world microscopic parameters affect 
the long-length-scale behavior. In this paper and in Q) 
we have provided oyery iew of the renormalization-group 
(|l67tll68tll69t llTOt llTl|) that many researchers use to un- 
derstand crackling noise and which is perhaps the most 
impressive use of abstraction in science. 

We now turn to some non-magnetic systems that have 
been particularly well studied. 



2. Martensites 

Hysteresis with crackling noise and in certain cases 
even microscopic return point memory has long been seen 
in fcrroclastic materials, such as shape memory alloys or 
martensites l|172[ Il73j) . A martensitic transformation is 
a diffusionless first order phase transformation where the 
lattice distortion is mainly described by a homogeneous 
shear. Many metals and alloys with a bec structure will 
upon cooling (or under strain) undergo this transition to 
a low temperature close packed structure. In athermal 
martensites, such as Cu-Zn-Al, thermal fluctuations are 
not relevant - temperature acts as an external driving 
field. Similar to ferromagnets ramping magnetic field, 
the martensitic transition takes place as a sequence of 
avalanches as the temperature is swept. The analogue 
of Barkhausen noise is actual noise: acoustic emission 



19 This subsection follow closely the presentation in reference y|), 
with some added references. 



20 These models are controversial Il38t ll.SHt) : they ignore catas- 
trophic external events like asteroids. 

21 The ra ndom -field model we describe was invented in this con- 
text l6 Ul40l) . as described in section ITT1 
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due to elastic waves in the ultrasonic range, generated 
by propagating domain walls during the avalanches. 

Experimental studies of the acoustic emission gener- 
ated during a thermally induced martensitic transfor- 
mation of a Cu-Zn-Al single crystal revealed an appar- 
ent absence of ch aracteristic scales in the distribution of 
avalanches 1J174T) . After an initial number of hysteresis 
cycles (ramping temperature) the system reaches a fi- 
nal attractor yielding two decades of power law scaling 
in the time distribution of acoustic emission pulse dura- 
tions. The authors argue that the cutoff at large times is 
probably due to experimental limitations (amplifier cut- 
off), and interpret the observed scaling behavior in terms 
of self-organized criticality. Of course, the system could 
have a tunable parameter that lies close to the critical 
point (two decades of power-law scaling arises in a range 
50% away from our critical point). 

In a more recent study, Perez-Reche, Stipcich, Vives, 
Manosa, Planes, and Morin (Jl75|) study the evolution of 
the kinetic features of the martensitic transition in a Cu- 
Al-Mn single crystal under thermal cycling. The authors 
use several experimental techniques, including optical mi- 
croscopy, calorimetry, and acoustic emission, to perform 
an analysis at multiple scales. Focussing especially on 
avalanches, they find that there are significant differences 
between the kinetics at large and small scales during the 
initial cycles. Upon repeated temperature cycles, how- 
ever, the system evolves from displaying a supercritical 
avalanche size distribution to a critical power-law dis- 
tribution. In the language of the RFIM this could be 
described as the effective disorder being increased by the 
temperature cycles from R < R c to R = R c . (The criti- 
cal exponents, and thus the universality class, is different 
from our RFIM). 

Both of these systems demand repeated cycling to pro- 
duce critical behavior. This is quite different from sys- 
tems with return-point memory, where repeated cycles 
don't change the behavior. A recently introduced, simple 
two dimensional phenomenological model for the marten- 
sitic transformation is shown to mimic near full reversal 
of morphology under thermal cycling in (J176U177T1 . 

A crossover from the hysteresis loops with a jump or 
"burst" to smooth hysteresi s loo ps (ramping tempera- 
ture) has been recorded in l|l78f ). as the grain size of 
macroscopic, polycrystalline specimen of an Fe-Ni-C al- 
loy is reduced. The crossover is explained in terms of fi- 
nite size effects in defect induced nucleation within each 
grain (due to extended defects like dislocation tangles, 
however, rather than random point disorder). It seems 
clear that one should look for critical fluctuations near 
the crossover that would promote an interpretation of the 
grain size as the analogue of the disorder parameter R in 
our model. 



3. Liquid Helium in Nuclepore 



In iJlJJl) Lilly, Finley, and Hallock show a hysteresis 
loop for the capillary condensation of superfluid helium 
in the porous material Nuclepore. Nuclepore is a polycar- 
bonate sheet perforated by a high density of nearly cylin- 
drical holes with a nominal diameter of about 2000A 
The pores intersect. The most probable number of in- 
tersections per pore is about 5. Most pores participate 
in a macroscopic interconnected network, (i.e. they are 
apparently internally connected beyond the percolation 
threshold). In the experiment the authors measure the 
amount of helium absorbed and capillary condensed in 
the Nuclepore sample as a function of the number of he- 
lium atoms present (chemical potential) at fixed temper- 
ature T = 1.52 K. 

Pore intersections result in weak interactions among 
the pores which lead to avalanches that have been ob- 
served upon slow withdrawal of helium from the cell. The 
upper branch of the hysteresis loop (upon withdrawal 
of He from the sample) actually has microscopic jumps. 
The largest avalanche observed involves 2.6 10 7 pores, 
which is about 2.6% of the system. The distribution of 
avalanche sizes is described by power laws. 

This system has clear analogies to our model. The cou- 
pling between the pores is certainly "ferromagnetic" : If a 
pore is emptied it becomes more likely for an intersecting 
pore to be emptied also. Different pore radii and different 
cross-sections of the pore intersections introduce a cer- 
tain randomness which might be used to tune the system 
to a critical point. In (180) the avalanches are studied 
for different characteristic pore sizes and the revelance 
of various theoretical model s wh ich one might apply to 
the system is discussed. In l)l8l|) a unifying theoretical 
approach to hysteresis phenomena associated with capil- 
lary condensation of gases in disodered mesoporous solids 
is given, and the connection with the zero-temperature 
theoretical studies of magnetization-reversal hysteresis of 
drive n in terfaces in disordered media is discussed (see 
also (182)). 



4. Superconductors 

a. Superconducting vortex avalanches In an interesting 
experiment on the dynamics of superconducting vortices 
in the Bean state in Nbir%Ti 53 %, Field, Witt and Nori 
recorded vortex avalanches as the system was driven to 
the threshold of instabi lity by the slow ramping of the 
external magnetic field (143). The individual avalanches 
were ranging in size from 50 vortices up to about 10 vor- 
tices, which is about 10 -6 of the total number of vortices 
in the system. The distribution of avalanches sizes was 
recorded for the events occurring in three 450 G wide 
windows of the external magnetic field range, centered 
at 2.25 kG, 5.33 kG, and 7.55 kG. At all three values for 
the field H there seems to be a power law up to some 
cutoff size, which varies with H. The authors suggest 
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that the self-organized criticality is the underlying mech- 
anism, despite the fact that their cutoff size changes as a 
function of H. 



b. Ultra-thin granular superconducting films in a parallel Held 
Wu and Adams discovered avalanches in film resistance 
corresponding to the collapse of macroscopic supercon- 
ducting regions involving up to 10 9 grains near the par- 
allel critical field in ultra-thin granular superconduct- 
ing Al- films l|183l) . The corresponding distribution of 
avalanche sizes displays four decades of power law scaling, 
(which in our simulations would correspond to a value 
of disorder about 10-15% away from the critical point). 
Fewer decades are seen at lower external magnetic fields 
— the field may act as a tunable parameter for a critical 
point. The authors point out the relation to our model 
and the zero temperature random bond Ising model with 
the same dynamics. The scaling exponent they extract 
from the power law decay is similar to the prediction 
from simulations of the two dimensional random bond 
Ising model. 



5. Prewetting on a disordered substrate 

Blossey, Kinoshita, and Dupont-Roc have studied the 
effect of substrate randomness on the dynamics of a 
prewetting transition l|184(K In their theoretical paper 
they first review several experiments and then suggest a 
mapping of the hysteretic prewetting system to a two di- 
mensional nonequilibrium RFIM in an external field H , 
by identifying the Ising spins ±1 with a wet or nonwet 
patch of the surface. They find that the covering of the 
substrate appears to be a critical transition. The growth 
exponents that they extract from their simulations are 
reported to be compatible with those of percolation, see 
also dl85h. 



X. UNSOLVED PROBLEMS 

There are many open and challenging questions in this 
field. 

1. Two Dimensions. Except for the experiments in 
magnets and martensites most of these experiments 
are in effectively two dimensional systems. D = 2 
appears to be the lower critical dimension of the 
transition in the RFIM Ij3t l4lf) . Extracting reliable 
scaling behavior in two dimensions remains a chal- 
lenge, and we have focused here on three dimen- 
sions where our understanding is reasonably com- 
plete. 

2. Avalanche average pulse shapes. The average 
pulse shapes predicted by the various theories (fig- 
ure !14|) do not appear close to those measured in ex- 
periments, which themselves disagree with one an- 
other ffigurc [TS|) . This is one of the few experiments 



where the scaling function has been extracted from 
the data. Are we mistaken about the universal- 
ity of these functions? Does this reflect measure- 
ment artifacts, or perhaps unrecognized important 
control parameters? Are there analytical correc- 
tions to scaling ( appendix IA"|) that we are missing 
that skew the curves? Or are there many dynami- 
cal universality classes? Recent experiments ( 186) 
suggest that the skewness of the smaller avalanches 
can be affected by external stress. Signatures of 
time asymmetries have also been seen i n high or - 
der spectra studies of Barkhausen Noise l)l87l : ll88f) . 



3. Correlated disorder and long range fields. 

Many of the experiments involve long-range fields 
which may in principle alter the universality class - 
perhaps in ways that differ from that of the infinite- 
range model. Some experiments involve different 
kinds of disorder, such as correlated disorder rather 
than point disorder and random anisotropics or 
random bonds rather than random fields (which 
probably will not change the universality class). 
Also, there are new possible critical behaviors when 
one starts with other local order parameters: for 
example, if we start not with Ising spins but with 
vectors (l9Sh . 



4. Extracting exponents and scaling functions. 

We need a more organized approach, both ex- 
perimentally and computationally, for extracting 
critical exponents and reporting scaling functions. 
These hysteresis models represent a serious chal- 
lenge, far more difficult than most of the critical 
points studied in equilibrium phase transitions. In 
part this is because the critical behavior is so much 
more important: because the avalanches become 
large so far from R c , one cannot get close to R c 
without having severe finite size effects - hence the 
R dependence of the parameters becomes a prob- 
lem. 
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APPENDIX A: Derivation of the scaling forms and 
corrections 

In our work we make extensive use of scaling col- 
lapses. 22 Many variations are important to us: Widom 
scaling, finite-size scaling, singular corrections to scaling, 
analytic corrections to scaling, rotating axes, and expo- 
nentially diverging correlation length scaling. The un- 
derlying theoretical framework for scaling is given by the 
renormalization group, in the context of equilibrium crit- 
ical phenomena and by now well explicated in a variety 
of texts. 

We have discovered that we can derive all the scaling 
forms and corrections that have been important to us 
from two simple hypotheses (found in critical regions): 
universality and invariance under reparameterizations. 
Universality is the statement that two completely differ- 
ent systems will behave the same near their critical point 
(for example, they can have exactly the same kinds of 
correlations) . Reparameterization invariance is the state- 
ment that smooth changes in the units or methods of 
measurement should not affect the critical properties - 
or rather, that the predictions of our scaling theories are 
only valid up to a smooth change of coordinates in the 
control parameters (the analytic corrections to scaling 
we discussed in section QJ . We use these properties to 
develop the scaling forms and corrections we use in this 
paper. Each example we cover will build on the previous 
ones while developing a new idea. 

Our first example repeats our discussion of size versus 
duration in section IIVI Consider some property F of a 
system at its critical point, as a function of a scale x. F 
might be the spin-spin correlation function as a function 
of distance x (or it might be the avalanche probability 
distribution as a function of size x, etc.) If two different 
experimental systems are at the same critical point, their 
.F's must agree. It would seem clear that they cannot be 
expected to be equal to one another: the overall scale 
of F and the scale of x will depend on the microscopic 
structure of the materials. The best one could imagine 
would be that 



Fi(a;i) = AF 2 {Bx 2 ) 



(Al) 



where A would give the ratio of, say, the squared mag- 
netic moment per domain of the two materials, and B 
gives the ratio of the domain sizes. 

Now, consider comparing a system with itself, but with 
a different measuring apparatus. Universality in this self- 
referential sense would imply F(x) — AF(Bx), for suit- 
able A and B. If instead of using finite constant A and 
B, we arrange for an infinitesimal change in the measure- 
ment of length scales, we find: 



F(x) = (1-ae) F[(l - e)x 



(A2) 



1 This appendix follows closely the presentation ir 



where e is small and a is some constant. Taking the 
derivative of both sides with respect to e and evaluating 
it at e = 0, we find —aF = xF', so 



F(x) ~x~ a . 



(A3) 



The function F is a power-law! The underlying reason 
why power-laws are seen at critical points is that power 
laws look the same at different scales. 

Now consider a new measurement with a distorted 



F 



(B(x) 



where 



measuring apparatus. Now F{x) ~ A 

A and B are some nonlinear functions. For example, 
one might measure the number of microscopic domains x 
flipped in an avalanche, or one might measure the total 
acoustic power B{x) emitted during the avalanche; these 
two "sizes" should roughly scale with one another, but 
nonlinear amplifications will occur while the spatial ex- 
tent of the avalanche is small compared to the wavelength 
of sound emitted: we expand B(x) — Bx + bo + b\/x + . . . 
Similarly, our microphone may be nonlinear at large 
sound amplitudes, or the absorption of sound in the 
medium may be nonlinear: A(F) = AF + a 2 F 2 + . . . 
So, 



A 



F 



(B(x)) 



A(F(Bx) + F'(Bx)(b + h/a 



F"{Bx) (...) + ••• + a 2 F z {Bx) 



(A4) 



We can certainly see that our assumption of universality 
cannot hold everywhere: for large F or small x the as- 
sumption of reparameterization invariance l|A4f) prevents 
any simple universal form. Where is universality pos- 
sible? We can take the power-law form F(x) ~ x~ a = 
x log A l log B which is the only form allowed by linear repa- 
rameterizations and plug it into (|A4|) , and we see that all 
these nonlinear corrections are subdominant (i.e., small) 
for large x and small F (presuming a > 0). If a > I , 
the leading correction is due to bo and we expect x~ a ~ x 
corrections to the universal power law at small distances; 
if < a < 1 the dominant correction is due to a 2 , and 
we expect corrections of order x~ 2a . Thus our assump- 
tions of universality and reparameterization invariance 
both lead us to the power-law scaling forms and inform 
us as to some expected deviations from these forms. No- 
tice that the simple rescaling led to the universal power- 
law predictions, and that the more complicated nonlinear 
rescalings taught us about the dominant corrections: this 
will keep happening with our other examples. 

For our second example, let us consider a property K 
of a system, as a function of some external parameter 
R, as we vary R through the critical point R c for the 
material (so r = R — R c is small). K might represent 
the second moment of the avalanche size distribution, 
where R would represent the value of the randomness; 
alternatively K might represent the fractional change in 
magnetization AM at the infinite avalanche ... If two 
different experimental systems are both near their critical 
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points (7*1 and r 2 both small), then universality demands 
that the dependence of K\ and K 2 on "temperature" R 
must agree, up to overall changes in scale. Thus, using a 
simple linear rescaling K(r) = (1 — fie)Ki (1 — e)r\ leads 
as above to the prediction 



K(r) = r 



(A5) 



Now let us consider nonlinear rescalings, somewhat dif- 
ferent than the one discussed above. In particular, the 
nonlincarity of our measurement of K can be dependent 
on r. So, 

Ar(K{r)\ = a a + a 1 r + a 2 r 2 + ... + a 01 K{r) + ... (A6) 

If \i > 0, these analytic corrections don't change the dom- 
inant power law near r — 0. However, if fi < 0, all the 
terms a n for n < —fj, will be more important than the sin- 
gular term! Only after fitting them to the data and sub- 
tracting them will the residual singularity be measurable. 
For the fractional change in magnetization: AM ~ r 13 
has < /3 < 1 (at least above three dimensions), so we 
might think we need to subtract off a constant term ao, 
but AM — for R > R c , so ao is zero. On the other 
hand, in a previous paper <|27|) . we discussed the singu- 
larity in the area of the hysteresis loop: Area ~ r 2 ~ a , 
where 2 — a = (3 + j38 is an analogue to the specific heat 
in thermal systems. Since a is near zero (slightly positive 
from our estimates of f3 and 5 in 3, 4, and 5 dimensions), 
measuring it would necessitate our fitting and subtract- 
ing three terms (constant, linear, and quadratic in r): we 
did not measure the area for that reason. 

For our third example, let's consider a function F(x, r), 
depending on both a scale x and an external parameter 
r. For example, F might be the probability Di nt that 
an avalanche of size x will occur during a hysteresis loop 
at disorder r = R — R c . Universality implies that two 
different systems must have the same F up to changes in 
scale, and therefore that F measured at one r must have 
the same form as if measured at a different r. To start 
with, we consider a simple linear rescaling: 



F 



(x, r) = (1 - ae) f((1 - e)x, (1 + Ce)r) . 



(A7) 



Taking the derivative of both sides with respect to e gives 
a partial differential equation that can be manipulated to 
show F has a scaling form. Instead, we change variables 
to a new variable y = x^r (which satisfies y' — y to order 
e). If F(x,y) = F(x,r) is our function measured in the 
new variables, then 



F{x, r) = F{x, y) = (1 - ae) F( (1 - e)x, y 



(A8) 



and -qF = xdF/dx shows that at fixed y, F ~ x~ a , 
with a coefficient J-(y) which can depend on y. Hence 
we get the scaling form 



This is just Widom scaling. The critical exponents a and 
C, and the scaling function J 7 (x l 'r) are universal (two 
different systems near their critical point will have the 
same critical exponents and scaling functions). We don't 
need to discuss corrections to scaling for this case, as 
they are similar to those discussed above and below (and 
because none were dominant in our cases). 

Notice that if we sit at the critical point r = 0, the 
above result reduces to equation (|A3|) so long as .F(O) is 
not zero or infinity. If, on the other hand, T(y) ~ y n 
as y —> 0, the two- variable scaling function gives a sin- 
gular correction to the power-law near the critical point: 
F(x,r) ~ x~ a T{x c r) ~ x" Q+ " c for x « r^ 1 ^: only 
when x ~ r~ 1 ''» will the power-law x~ a be observed. 
This is what happened in two dimensions to the in- 
tegrated avalanche size distribution and the avalanche 
correlation functions (figures 35, 36, and 37b in refer- 
ence I©). 

For the fourth example, we address finite-size scaling of 
a property K of the system, as we vary a parameter r. If 
we measure K(r, L) for a variety of sizes L (say, all with 
periodic boundary conditions), we expect (in complete 
analogy to (|X9l) ) 



K^^-r-t 1 K,{rL l l u ). 



(A10) 



F(x,r) 



T(x 



C 



(A9) 



Now, suppose our "thermometer" measuring r is weakly 
size-dependent, so the measured variable is C(r) = r + 
c/L + C2/L 2 + . . . The effects on the scaling function is 

K(c(r),L\ ~ r ->* x 

(KirL 1 '") + 
(cL 1 ^- 1 + c 2 L 1/u - 2 ) JC'(rL 1/l/ ) + ...). (All) 

In two and three dimensions, v > 1 and these correc- 
tion terms are subdominant. In four and five dimensions, 
we find 1/2 < v < 1, so we should include the term mul- 
tiplied by c in equation (|A11|I . However, we believe this 
first term is zero for our problem. For a fixed boundary 
problem (all spins "up" at the boundary) with a first or- 
der transition, there is indeed a term like c/L in r(L). 
At a critical transition, the leading correction to r(L) 
can be c/L or a higher power of L {1/L 2 and so on). 
This seems to depend on the model studied, the geom- 
etry of the system, and the boundary conditions (free, 
periodic, ferromagnetic, . . .) Furthermore, for the same 
kind of model, the coefficient c itself depends on the ge- 
ometry and boundary conditions, and it can even vanish, 
which leaves only higher order corrections. In a periodic 
boundary conditions problem like ours, we expect that 
the correction is smaller than c/L. Our finite-size scaling 
collapses for spanning avalanches N, the second moments 
(S 2 ), and the magnetization jump AM, were successfully 
done by letting c = 0. 

For the fifth example, consider a property K depending 
on two external parameters: r (the disorder for example) 



2S 



and h (could be the external magnetic field H 
Analogous to (|A9() . K should then scale as 



K{r,h) r^r-i* ICih/r 136 ) 



H r 



(A12) 



Consider now the likely dependence of the field h on the 
disorder r. A typical system will have a measured field 
which depends on the randomness: C(h) = h+br+b2r 2 + 
. . . (Corresponding nonlincaritics in the effective value of 
r are subdominant.) This system will have 



K 



(r,C(h)\ = r-i* x 

(lC{h/r fiS ) + {br + b 2 r 2 ) r~ ps /C'(/i/^ 5 ))(A13) 



Now, for our system 1 < [35 < 2 for dimensions three and 
above. This means that the term multiplied by b is domi- 
nant over the critical scaling singularity: unless one shifts 
the measured h to the appropriate hi = h + br, the curves 
will not collapse (e.g., the peaks will not line up horizon- 
tally). We measure this (non-universal) constant for our 
system using the derivative of the magnetization with 
field dM/dH(r,h). The magnetization M(r,h) and the 
correlation length £(r, h) should also collapse according 
to equation (|A12|1 (but with h+br instead of ft); we don't 
directly measure the correlation length, and the collapse 
of M(r, h) in the bottom figure IT51 includes the effects of 
the tilt b. In two dimensions, (35 is large (probably infi- 
nite), so in principle we should need an infinite number of 
correction terms: in practice, we tried lining up the peaks 
in the curves (with no correction terms); because we did 
not know (3 (which we usually obtained from AM, which 
gives (3jv = in two dimensions), we failed to extract 
reliable exponents in two dimensions from dM/dH. 

For the sixth example, suppose F depends on r, h, and 
a size x. Then from the previous analysis, we expect 



F{x,r,h) ~x~ a f(iV, h/r p ") 



(A14) 



Notice that universality only removes one variable from 
the scaling form. One could in practice do two-variable 
scaling collapses (and we believe someone has probably 
done it), but for our purposes these more general scaling 
forms are used by fixing one of the variables. For exam- 
ple, we measure the avalanche size distribution at various 
values of h (binned in small ranges), at the critical dis- 
order r = 0. We can make sense of equation <|A14(1 by 
changing variables from h/r^ s to x**P 5 h: 



F(x,r,h) ~x~ a f(x c r, x^ s h). 



(A15) 



Before we can set r — 0, we must see what are the pos- 
sible corrections to scaling in this case. If the disorder r 
depends on the field, then instead of the variable r, we 
must use r + ah (the analysis is analogous to the one in 
example five; other corrections are subdominant). Set- 
ting r = now, leaves F dependent on its first variable, 
as well as the second: 

F(x,r,h) ~ x~ a f(x c (ah), x c0s h) w x~ a x 



(#(0, xfWh) + 

ahx< T^°H0, x^ 5 h)), (A16) 



where ft 1 ' ) is the derivative of T with respect to the 
first variable (keeping the second fixed). 

For the binned avalanche size distribution, x^ is S a , 
where 0<er<l/2as we move from two dimensions to 
five and above. Thus, the correction term will only be 
important for rather large avalanches, S > h~ 1 ' cr , so long 
as we are close to the critical point. Expressed in terms 
of the scaling variable, important corrections to scaling 
occur if the scaling variable X = S a/3S h > h 1 " 136 . For us, 
(36 > 3/2, and we only use fields near the critical field 
(h < 0.08), so the corrections will become of order one 
when X = 4 for the largest h we use. In 3 and 4 di- 
mensions, this correction does not affect our scaling col- 
lapses, while in 5 dimensions some of the data needs this 
correction. We have tried to avoid this problem (since we 
don't measure our data such that it can be used in a two- 
variable scaling collapse) by concentrating on collapsing 
the regions in our data curves where this correction is 
negligible. 

A similar analysis can be done for the avalanche time 
distribution, which has two "sizes" S and t and one pa- 
rameter r which is set to zero; because we integrate over 
the field h the correction in l|A16|) does not occur, and 
other scaling corrections are small. 

Finally, we discuss the unusual exponential scaling 
forms we developed to collapse our data in two dimen- 
sions. If we assume that the critical disorder R c is zero 
and that the linear term in the rescaling of r vanishes 
(£er in equation (|A7() vanishes), then from symmetry the 
correction has to be cubic, and equation (|A7|I becomes: 

F(x,r) = (1-ae) F^l-e)^, [l + ker 2 )rY (A17) 

with k (which is not universal) and a constants, and e 
small. 

Taking the derivative of both sides with respect to 
e and setting it equal to zero gives a partial differen- 
tial equation for the function F. To solve for F, we 
do a change of variable: (x, r) — > (x, y) with y = 
x e~ a ' r . The constant a* is determined by requiring 
that y rescales onto itself to order e: we find a* = 1/2 k. 
We then have: 



dF 

-a F(x,y) -— x 



which gives 



F{x,r)=x- a f(xe- 1/2kr2 y 



(A18) 



(A19) 



This is one of the forms we use in 2 dimensions for 
the scaling collapse of the second moments (S 2 ) int , the 
avalanche size distribution Di nt integrated over the field 
H, the avalanche correlation Gi„t, and the spanning 



2!) 



avalanches N. We use another form too which is ob- 
tained by assuming that the critical disorder R c is not 
zero but that the linear term in the rescaling of r still 
vanishes. Instead of equation l|A17(l . we have: 

F(x,r) = (l-cee) F((l-e)x, (l + ier)r). (A20) 

The function F becomes: 

F(x, r) = x- a T(xe-^ £r ) . (A21) 

The corrections to scaling for the last two forms (equa- 
tions l|A19|l and l|A21jl ) are similar to the ones discussed 
above. They are all are subdominant. 
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